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Program  activities  tc  accomplish  these  uLjectives  were  organized  "nder  five 
tasks.  Under  Task  1  -  Literature  Survey,  a  survey  was  performed  to  dfci.'rmlne 
the  state-of-the-art  In  design  .and  analysis  of  bolted  composite  Joints. 
Experimental  evaluations  of  joint  static  strength  were  performed  under  Tasks  2 
and  3.  In  Task  2  -  Evaluation  of  Joint  Design  Variables,  strength  data  were 
obtained  through  an  experimental  program  to  evaluate  the  effects  of  twelve 
joint  design  variables.  In  Task  3  -  Evaluation  of  Manufacturing  and  Service. 
Anomalies^  effects  of  seven  anomalies  on  joint  strength  were  evaluated  experi¬ 
mentally  and  compared  with  Task  2  strength  data.  Bolted  composite  joint 
durability  was  evaluated  under  Task  4  -  Evaluation  of  Critical  Joint  Design 
Variables  on  Fatigue  Life.  Seven  critical  design  variables  or  manufacturing 
anomalies  were  Identified  based  on  Task  2  and  3  atier.;''.).'  data.  Under  Task  5  - 
Final  Analyses  and  Correlation,  required  data  reduction,  methodology  development 
and  correlation,  .and  necessary  documentation  were  performed. 

This  report  documents  all  program  activ.itle.s  pertoi’med  under  Tasks  2,  3, 

4  and  5.  Activities  performed  under  Task  1  -  Liter.. ture  Survey,  were  pre¬ 
viously  reported  on  AFFDL-TR-78-179. ^Static  strength  methodology  and  evaluations 
of  joint  static  and  fatigue  test  data  are  reported.  Analytic  studies  complement 
mcithodology  development  and  Illustrate:  the  need  for  detailed  stress  analysis, 

I  the  utility  of  the  developed  "Bolted  Joint  Stress  Field  Model"  (BJSFM)  procedure, 
and  define  model  limitations.  For  static  strength  data,  correlations  with 
analytic  predictions  are  Included.  Data  trends  In  all  cases  are  discussed 
relative  to  joint  strength  and  failure  node.  For  joint  fatigue  studies,  data 
trends  are  discussed  relative  to  life,  hole  elongation,  and  failure  mode 
behavior. 

This  final  report  is  organized  Ju  the  folxowlng  three  volumes: 

/  Volume  1  -  Methodology  Development  and  Data  Evaluation 

Volume  2  Test  Data,  Equipment  and  Procedures 

Volume  3  -  Bolted  Joint  Stress  Field  Model  (BJSFM)  Computer  Program 
User ' 8  Manual 
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SECTION  I 


INTRODUCTION 


One  objective  of  this  program  was  to  develop  a  static 
strength  methodology  for  mechanically- fastened  composite  joints. 
Tills  volume  documents  user-options  and  instructions  for  a  com¬ 
puter  program  to  analyz€>  the  effects  of  stress  concentrations  on 
laminate  strength.  Entitled  "Bolted  Joint  Stress  Field  Model" 
(BJSPM',  it  computes  stress  distributions  on  a  lamina  or  laminate 
basis  for  unloaded  or  loaded  (bolt  bearing)  holes  in  isotropic  or 
anisotropic  materials.  Failure  predictions  based  on  leunlna 
properties  and  one  of  several  failure  criteria  are  possible. 
This  volume  describes  the  formulation  and  input  data  requirements 
and  output  options .  Sample  problems  and  a  computer  program 
listing  are  also  included. 


SECTION  II 


PROGRAM  DESCRIPTION 


The  Bolted  Joint  Stress  F’ield  Model  has  been  developed  to 
facilitate  strength  analysis  of  isotropic  or  anisotropic 
materials  at  individual  fastenei'  holes.  Static  str*^nqth  of  an 
anisotropic  laminate  with  a  fastener  liole  is  predicted  using  a 
closed- form  analytic  approach  based  on  (1)  elastic  anisotropic 
theory  of  elasticity,  (2)  lamination  plate  theory  and  (3)  one  of 
several  optional  failure  hypotheses.  The  program  has  capability 
to  handle  strength  and  stiffness  anisotropy,  general  in-plane 
loadings,  as  shown  in  Figure  1,  multi-material  (hybrid)  laminates 
and  arbitrary  hole  (bolt)  sizes.  BJSFM  modular  substructuring  is 
illustrated  in  Figure  2.  Input  data  required  are:  lamina  mech¬ 
anical  properties,  in-plane  loadings,  hole  geometry,  and  hole 
loading.  Options  are  available  which  provide  computation  results 
after  each  program  block. 

The  stress  field  calculations  are  based  on  two-dimensional 
anisotropic  theory  of  elasticity  solutions  for  a  homogeneous, 
anisotropic  infinite  plate.  Laminate  stress  distributions  arovmd 
an  unloaded  or  loaned  (bolt  bearing)  hole  are  calculated  using 
plane  stress  assumptions.  Laminate  stress  and  strain  di8trib\i- 
tions  for  combined  bearing  and  bypass  loads  are  obtained  using 
the  principle  of  superposition  (Figure  3).  Fastener  bearing  in 
idealized  as  a  cosine  radial  stress  distribution  (Figure  4). 
Finite  width  corrections  for  loaded  holes  are  based  on  super¬ 
position  of  infinite  plate  results  as  shown  in  Figure  5. 
Infinite  plate  solutions  are  exact  while  coirections  for  finite 
width  joints  are  approximate  and  most  accurate  for  width- to- 
diameter  ratios  greater  than  four. 

Iiaminate  strains  are  calculated  using  material  compliance 
relations.  Laminate  compliance  coefficients  are  determined  using 
lamination  plate  theory  with  unidirectional  (lamina)  elastic 
constants,  lamina  orientations  and  thicknesses.  Strains  for 
individual  plies  along  lamina  principal  material  axes  are  calcu¬ 
lated  using  coordinate  transformations.  The  solution  is  strictly 
valid  only  for  homogeneous  media;  however,  it  has  been  assumed 
valid  for  mid-plane  symmetric  laminates. 

Laminate  failure  is  predicted  by  comparing  elastic  stress 
distributions  with  any  of  five  material  failure  criteria  on  a 
ply-by-ply  basis.  Failure  can  be  assessed  at  any  location  in  the 
field  of  the  plate. 
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Figur«  1.  Q«n«raS  Load  Conditiona  Anaiyzad  UsinQ  SJSFm 


FIgura  2.  Boltad  Joint  Straaa  Flold  Modal  -  BJSFM 

Computer  Program  Flow  Chart 
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FIgur*  3.  Superposition  of  Linoar-EI««tlc  Stross  Solutions 
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FIgurs  4.  Assumod  Cosins  Bolt-Losd  Distribution 
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Figure  5.  Superposition  of  Infinite  Piste  Results  Approximates  a  Finite  Width  Joint 
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SECTION  III 


USER'S  INSTRUCTIONS 


Conversational  interactive  procedures  are  used  for  specify¬ 
ing  input  data  for  operation  of  the  BJSPM  program  which  is  pro¬ 
grammed  to  accept  free  formatted  inpuu  data.  A  user  may,  after 
becoming  familiar  with  the  input  procedures,  elect  to  delete 
input  instructions  and  receive  only  question  marks,  identifying 
all  required  input  data.  Using  the  various  output  options,  users 
may  receive  as  much  or  as  little  data  as  desired.  The  nature  of 
required  input  data  is  dependent  on  user-selected  output  data 
options;  the  BJSFM  program  automatically  adjusts  its  input  data 
requirements  to  accommodate  each  output  data  option. 

Up  to  eight  different  ply  orientations  and  three  different 
materials  may  be  input.  For  each  ply  orientation,  a  correspond¬ 
ing  thickness  must  be  specified  as  well  as  material  for  hybrid 
laminates.  Ply  thickness  may  be  either  actual  or  a  normalized 
thickness.  A  mid-plane  symmetric  stacking  sequence  is  assumed. 
Zero  degree  plies  are  oriented  parallel  to  the  X-axis.  Nomencla¬ 
ture  applicable  to  tl^e  BJSFM  is  summarized  in  Figure  6.  Positive 
angles  are  measured  Counterclockwise  from  the  X-axis.  All  input 
data  units  must  be  consistent. 

Unidirectional  lamina  material  stiffness  properties  aro 
required  input  data  for  each  different  material  specified. 
Unidirectional  lamina  strength  allowables  for  each  material  are 
required  only  if  failure  analysis  is  to  be  performed.  If  the 
maximum  strain  material  failure  criterion  is  used,  lamina  strain 
allowables  must  be  input;  otherwise,  input  leunina  allowables  are 
in  terms  of  stress. 

Any  set  of  in-plane  far- field  stresses  may  be  applied  to  an 
infinite  anisotropic  or  isotropic  plate  (Figure  1).  Bearing 
stress  direction  is  independent  of  far-field  stress  directions. 

The  BJSFM  is  only  capable  of  handling  finite  widths  for  bolt 
bearing  problems;  the  width,  W,  is  defined  as  perpendicular  to 
the  bolt  load  direction  (Figure  1).  The  stresses  calculated  in  a 
finite-width  bolt  bearing  problem  are  approximate  and  most 
accurate  for  width-to-dieimeter  ratios  greater  than  four.  In 
combined  loading  conditions,  the  finite  width  routine  applies  to 
only  the  loaded  hole  portion  of  the  problem.  To  obtain  infinite 
plate  results  for  a  loaded  hole,  input  specimen  width  as  0.0. 

The  user  must  specify  the  "range"  (between  low  and  high)  of 
angular  interval  between  locations  around  the  hole  for  which  data 
will  be  calculated.  This  range  must  also  be  su.xlivided  by  user- 
selected  "degrees  between  output"  to  specify  points  at  which 
calculations  are  to  be  made.  "Step  increments"  are  used  to 
obtain  data  at  increasing  distances  away  from  the  hole  boundary. 
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Coordilwt*  Systanw 

X  -  Y  Laminate  Axis  System  Originating  at  the  Center  of  the  Hole.  Zero  Degree  Plies  are  Parallel  to  the 
X-Axis. 

X'  -  Y'  Rotation  of  the  X-Y  Axis  System  for  Application  of  Far-Field  Stresses 
1  -  2  Lamina  Axis  System;  Fibers  are  Parallel  to  the  1  Axis  and  Transverse  to  the  2  Axis. 

Variable  Dascription  (Units)** 

El  Lamina  Modulus  of  Elasticity  in  Fiber  (1)  Direction  (F/L  ) 

E2  Lamina  Modulus  of  Elasticity  in  Transverse  (2)  Direction  (F/L  ) 

G12  Lamina  Shear  Modulus  (F/L^l 

V12  Lamina  Poisson's  Ratio 

EX  Laminate  Modulus  of  Elasticity  in  X  Direction  (F/L  ) 

EY  Laminate  Modulus  of  Elasticity  in  Y  Direction  (F/L^) 

SXY  (.aminate  Shear  Modulus  (F/L^l 

VXY  Laminate  Poisson's  Ratio 

T1  Lamina  Allowable  Tensile  Strength  in  Fiber  (1)  Direction  (F/L  or  L/L) 

Cl  Lamina  Allowable  Compressive  Strength  in  Fiber  (1)  Direction  (F/L^  or  L/L) 

T2  Lamina  Allowable  Tensile  Strength  in  Transverse  (2)  Direction  (F/L^  or  L/L) 

C2  Lamina  Allowable  Compressive  Strength  in  Transverse  (2)  Direction  (F/L^  or  L/L) 

S  l  amina  Allowable  Shear  Strength  (F/L^  or  L/L) 

PX  Stress  in  X'  Direction  (F/L^)  -  Independent  of  Input  Thickness 

PY  Stress  in  Y'  Direction  (F/L^)  •  Independent  of  Input  Thickness 

PXY  Shear  Stress  (F/L^)  •  Independent  of  Input  Thickness 
P  Applied  Bearing  Stress  (F/L^!  •  P  =  Bolt  Load/(Dia  x  Actual  Thickness) 

U  Displacement  in  X  Direction  (L) 

V  Displacement  in  Y  Direction  (L) 

W  Specimen  Width  (L)  -  Bolt  Loading  Only 

DIST  Radial  Distance  from  Hole  Boundary  (L) 

a  Aiigle  of  Applied  Bolt  Load  with  X  Axis  (i.e..  Bolt  Loading  Angle) 

P  Rotation  Angie  of  X'  •  Y'  Axes  from  X  -  Y  Axis  System  (i.e.,  Off  Axis  Loading  Angle) 

0  Angle  from  X  Axis  to  a  Point  Around  Fastener  Hole 

0  Rotation  Angle  of  1  -  2  Axes  from  X  •  Y  Axis  System  (i.e..  Ply  Orientation  Angle) 

All  angular  meaiuremantt  are  positive  counterclockwise  from  the  X  axis. 

*  See  also  Figure  1 . 

**  Any  consistent  sat  of  units  may  be  used: 

F  •  Force 
L  ■  Length 

aFi»aiis-22S 

Figurv  6.  Program  Nomenclature* 
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The  maximum  input  step  increment  io  seven  evenly  spaced  concen¬ 
tric  circles;  the  first  step  is  always  at  the  hole  boundary. 


The  option  to  use  any  one  of  five  different  failure  criteria 
has  been  programmed  into  the  BJSFM.  Failure  analysis  is  applied 
on  a  ply-by-ply  basis;  therefore,  only  unidirectional  (lamina) 
allowables  are  required  input  data.  Only  the  maximum  strain  cri¬ 
terion  requires  the  allowables  to  be  input  as  strains;  all  others 
jse  stress  allowables.  Equations  for  each  of  the  programmed 
failure  criteria  are  given  below.  When  the  right  hand  side  of 
any  of  the  equations  exceeds  unity,  failure  has  been  predicted 
for  the  ply.  Tension  or  compression  stress/strain  allowables 
used  in  each  criteria  are  selected  automatically,  depending  on 
the  sign  of  individual  stress  field  components  being  evaluated. 


Maximum  Strain 


1 


Maximum  Stress 


1 


Tsai-Hill 


Modified  Tsai-Wu 


Hoffman 


An  example  printout  of  the  “conversational" 
request  input  data  is  shown  in  Figure  7. 


language  used  to 


RHH 

IDO  YOU  WANT  INSTRUCTIONS? 

PYES 

SELECT  DESIRED  OUTPUT  FROfI  THE  FOLLOV/ING  CASES. 

1  CARPET  PLOT  DATA 

2  LAMINATE  PROPERTIES 

3  LAMINATE  STRESSES 

4  LAMINATE  STRAINS 

5  CIRCUMFERENTIAL  Ci  RADIAL  STRESSES/ STRAINS 

6  DISPLACEMENTS 

7  STRAINS  PER  PLY 
9  STRESSES  PER  PLY 

9  FAILURE  CRITERIA  PER  PLY 
10  AUTOMATIC  SEARCH  FOR  FAILURE 
72,3,4,5,6,7,8,9,10 

INPUT  NUMBER  OF  DIFFERENT  PLIES  TO  BE  INPUT  (8  MAX)  AND 
NUMBER  OF  DIFFERENT  MATERIALS  (3  MAX) 

74,  1 

INPUT  THE  UNIDIRECTIONAL  MATERIAL  PROPERTIES  FOR  EACH  MATERIAL 
IN  THE  FOLLOWING  ORDER:  El,  E2,  G12,  POISSONS  RATIO 
718.8586, 1.9E6, .85E6, .3 

INPUT  THE  UNIDIRECTIONAL  ALrXR7ABLES  FOR  EACH  MATERIAL 
IN  THE  FOLIXXaNG  ORDER:  Tl,  Cl,  T2,  C2,  SHEAR 
7230000, 320000, 20200,32300, 17300 
INPUT  THE  ANGULAR  ORIENTATION  OF  EACH  PLY 

70. . 45. ,-45. ,90. 

INPUT  THE  THICKNESS  OF  EAOi  PLY 

7.5. . 2. .2. .1 

INPUT:  FAR  FIELD  STRESSES  PX,PY,PXY,  OFF  \XIS  ANGLE,  BEARING  S\  TSS 

AND  DOLT  LOADING  ANGLE. 

710000 . .  0 . . 2500 . . 45 . . 50000 . . 10 . 

I'TpuT  WIDTH  (0.0  FOR  INFINITE  PLATE) 

71.5 

INPUT  BOLT  DIAMETER,  DEGREES  BETWEEN  OUTPUT,  LO»/  RANGE.  HIGH  RANGE, 
STEP  INCRDIENT  AND  NUrlDER  OF  STEPS  DESIRED  (7  MAX) 

7. 25. 30. . 0.. 360.. .02.2 

INPUT  THE  NUMBER  V/HICH  CORRESPONDS  TO  THE  FAILURE  CRITERIA 
YOU  V/ISH  TO  USE 

1  MAXIMUM  STRAI N 

2  MAXIMUM  STRESS 

3  TSAI -HILL 

4  MODIFIED  TSAI-vn; 

5  HOFFMAN 
73 


Figurt  7.  Exampls  input  Instructions 
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SECTION  IV 


OUTPUT  OPTIONS 


Various  output  options  are  available  for  user  selection. 
The  user  may  select  any  or  all  of  the  following  options  by  input¬ 
ting  the  appropriate  number(s): 

1.  Carpet  Plot  Data 

2 .  Laminate  Properties 

3.  Laminate  Stresses 

4 .  Laminate  Strains 

3.  Circumferential  and  Radial  Stresses/Strains 

6.  Laminate  Displacements 

7.  Strains  per  Ply 

8.  Stresses  per  Ply 

9.  Failure  Criteria  per  Ply 

10.  Automatic  Search  for  Failure 

A  brief  description  of  each  output  option  follows. 

Option  #1  -  Carpet  Plot  Data  -  The  carpet  plot  data  routine 
will  automatically  vary  the  layup  of  a  user  input  0*/+i}i/90'’ 
leuninate  feunily  and  calculate  any  one  or  ell  of  the  other  output 
options  (2  through  10)  for  each  layup.  Sixty-six  different 
layups  are  automatically  calculated  in  thit  routine;  therefore, 
large  amounts  of  data  will  be  generated  when  using  this  output 
option. 

Option  #2  -  Laminate  Properties  -  Laminate  stiffness  proper¬ 
ties  are  calculated  using  the  unidirectional  material  elastic 
constants,  ply  angular  orientations  and  ply  thicknesses.  These 
properties  are  calculated  with  respect  to  the  X-Y  axes  and  are 
the  same  as  would  be  obtained  using  conventional  lamination 
theory  approaches . 

Options  #3  and  #4  -  Laminate  Stresses/Strains  -  Laminate 

stress  and  strain  distributions  are  available  as  output  at  points 
around  the  perimeter  of  the  hole  and  at  other  user-specified 
concentric  circles  about  the  hole  boundary.  Principal  stresses 
and  strains  are  also  calculated.  All  output  is  referenced  to  the 
X-Y  axes.  Points  are  locatesd  by  the  radial  distance  away  from 
the  hole  boundary  and  the  angular  orientation  from  the  X-axis 
(Figure  8) . 

Option  #5  -  Circumferential  and  Radial  Stresses/Strains  - 
Circumferential  and  radial  laminate  stresses  and  strains  are 
calculated  by  a  coordinate  transformation.  Output  is  in  polar 
coordinates . 
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LAMINATE  STRESSES 


DI3T 

ANGI£ 

X  STKESS 

Y  STRESS 

0.000 

0.00 

-62603.22 

44412.49 

0.000 

30.00 

-40151.39 

-on. 39 

0 . 000 

60.00 

-6019.06 

-29547.53 

0 . 000 

90.00 

50351 .00 

-11003.21 

0.000 

120.00 

401 27.47 

16033.31 

L 

150.00 

9017.46 

70047.34 

■■■  — _ 

-  t.*-"  1  ^ 

SHEAR 

STRESS 

-9.44 

-34070.37 

-19090.77 

-9.44 

77703.47 


MAX. 


HIN. 


PRINCIPAL  PRINCIPAL 


44417.50 
10054.04 
4557.17 
5035* 


-67093.73 

-59074.32 


DIRECTION 
.01 


Figure  8.  Laminate  Street  Dietribution  Output  Data  Option 


Option  #6  -  Laminate  Diaplacements  -  Displacements  for  each 
point  are  output  as  U  and  V,  which  are  displacements  in  the  X  and 
Y  directions  reEpectively .  Due  to  limitations  in  the  derivation, 
displacements  for  the  loaded  hole  case  shall  be  considered  accu¬ 
rate  only  within  approximately  three  times  the  diameter  of  the 
fastener .  Unloaded  hole  displacements  are  exact  throughout  the 
plate . 

Options  #7  and  8  -  Strains/stresses  per  Ply  -  Strains  and 
stresses  per  ply  are  calculated  and  output  in  the  launina  (1-2) 
coordinate  system.  Eacli  ply  is  identified  along  with  the  loca¬ 
tion  of  the  point  around  the  hole  for  which  stresses/strains  ar* 
calculated  (Figure  9). 
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STRAINS  PER  PLY 


DIST 

ANGLE 

PLY 

STRAIN  1 

STRAIN  2 

SHE^AR  STRA 

0.000 

0.00 

0.00 

-.007128 

.011127 

— . 000004 

0.000 

0.00 

A5.00 

.001990 

.002001 

.010254 

0.000 

0.00 

-4^.00 

.002001 

.001990 

-.0107^ 

0.000 

0.00 

90 . 00 

.011127 

_ 

0.000 

30.00 

0.00 

-.003509 

OPIMIIMM 


Flgur*  8.  Pljf'by'Ply  Strain  Distribution  Output  Data  Option 


Option  #9  -  Failure  Criterici  per  Ply  -  Failure  criteria  per 
ply  option  applies  the  user-selected  material  failure  criterion 
(max.  stress,  Tsai-Hill,  etc.)  to  each  ply  using  the  input 
material  allowables.  The  "failure  number"  obtained  as  output 
data  indicates  the  value  calcvilated  by  the  failure  criterion 
using  the  stress  or  strain  components  at  a  point.  A  failure 
number  equal  to  or  greater  than  one  predicts  ply  failure.  The 
progreun  automatically  selects  tension  or  compression  allowables 
depending  on  the  sign  of  individual  sUress/strain  components 
being  evaluated.  "Failure  ratios"  are  output  v/hich  indicate  the 
relative  magnitude  of  contributing  stress  ciwiponents  to  the 
overall  failure  number.  Therefore,  failure  can  be  assessed  as  to 
which  stress  component  is  most  signifies  it.  These  failure  ratios 
are  in  terms  of  the  Icimina  (1-2)  coordinate  system. 
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Option  #.10  -  Aut'Omatic  Search  for  Failure  -  The  automatic 
searcn  For  FaTIui e  routine  will  uearch  over  a  uaer-epecif ied 
range  aV.  each  am  ular  inclement  for  the  moet  critical  single 
p>oint  as  calculated  by  the  material  failure  criterion.  Search 
for  failure  is  onlv  done  at  the  first  step  increment  away  from 
the  hvole  boundary.  Therefore,  if  a  search  for  failure  is  to  be 
performed  at  the  bouAdary  of  the  hole,  the  step  increment  must  be 
input  as  0.0.  The  program  will  automatically  ratio  the  input 
atreus  field  until  first  ply  failure  is  predicted.  Output  is  the 
in-plane  stresses  at  vhich  failure  is  predicted  along  with  the 
angular  orientation  ol'  the  predicted  failure  location.  Failure 
numbers  are  also  given  for  all  other  plies  at  thva  critical  ply 
failure  initiation  angle.  Failure  ratios  are  also  output. 
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SECTXON  V 


PROGRAM  LIMITATIONS 


The  following  are  the  limitations  of  the  BJSFM  program. 

o  Strictly  valid  for  homogensoue  anisotropic  flat  plates 
and  assumed  valid  for  mid-plane  synmetric  laminates. 

o  Displacements  inaccurate  at  points  more  than  three  times 
the  hole  diameter  away  from  the  hole  boundary  for  loaded 
hole  cases. 

o  Stress  fields  inaccurate  for  width  -  Lo-diameter  ratios 
less  than  four. 

o  Maximum  of  eight  different  ply  angular  orientations 
( input) . 

o  Maximum  of  three  different  materials  for  hybrid  laminates 
{ input) . 

o  Maximum  of  seven  steps  away  from  the  hole  (output) . 


The  following  equation  must  be  satisfied  to  obtain  valid  output. 
[(High  Range) “(Low  Range)] /(Degrees  Between  Output)  £72 


The  following  data  must  be  input  as  integers: 

Output  Option  Numbers 
Number  of  Different  Plies 
Number  of  Different  Materials 
Material  Number 
Number  of  Steps 
Failure  Criteria  Number 


loo  YOU  WANT  IMtTRUCTIOMS? 

TTM 

sa.Krr  ocaired  outpitt  prom  thi  polumimo  cams. 

l  CARPET  PLOT  DATA 
LAHIMATB  PROPCRTII8 
LAMINATE  STRESSES 
LAMINATE  STRAINS 

CXRCUMPEREHTIAL  4  RATIAL  STRESSES/ STRAINS 
OISPLACENBtlTS 
STRAINS  PER  PLY 
STRESSES  PER  PLY 
PAIL’mr  CRITERIA  PER  PLY 
AUTOMATIC  SEARCH  POR  PAtrAlRS 


I 
I 

4 

5 

s 

7 
t 

» 

10 
74, S 

INPUT  NUMBER  OP  OIPPERBHT  PLIES  TO  BE  INPUT  (S  MAX)  AND 
NUMBER  or  DirrCRENT  MATERIALS  (1  MAX) 

- M,l-| 

liWr  THE  UNIDIRECTIONAL  MATERIAL  PROPERTIES  POR  SAC  MATERIAL 
IN  THE  rOLLOUIlia  OROCRi  Cl,  BS,  01},  POISSONS  RATIO 
*— nr,. R5C6,1.<»C«,  .8SS4.  .1 

INPUT  iMC  ANGULAR  ORIENTATION  OP  BACH  PLY 
■70. .4S., -44. ,90. 


r—  70. ,95 
— j  INPUT 

•--7.42,  . 


PLY 


THE  THICKNESS  OP  BACH 
2S..25,.00  - 

IHPUTi  PAR  PIELD  STRESSES  PX,PY,PXY. 


1-  42%  -  o"  Plies,  50%  ±45°,  8%  90° 


OPP  AXIS  ANOLB,  BBARINO  STRBSS 

AND  BOLT  LOADIMO  ANOLB. 

70,0,0,0,40000. ,10. 

INPUT  WIDTH  (0.0  POR  INPIHITB  PLATE) 

71.4 

INPUT  BOLT  DIAMETER,  DBORCBS  BETWEEN  OUTPUT,  LOW  RANOB,  HIGH  RANOE, 
STEP  IMCRMENT  AND  NUMBER  OP  STEPS  DBSIRED  (7  MAX) 

7. 24, 10. ,0., 90., .02,3 


r  ' — 

LAMINATE 

STRAINS  (Output  Option  No.  4) 

DIST 

ANGLE 

X  STRAIN 

Y  STRAIN 

SHEAR 

MAX. 

Hit). 

DIRBCTIC 

STRAIN 

PRt^K:IPAL 

PRINCIPAL 

o.roo 

-  0.00 

-.006927 

.007109 

. onoool 

.007109 

-.006977 

-.00 

0.000 

10.00  1 

-.004117 

.001918 

-.017140 

.005704 

-.008004 

18.20 

0.000 

Lrso.oo 

.001704 

-.006481 

-.oloiin 

.004099 

-.008979 

-12.87 

0.000 

•-90.00-' 

.004412 

-.000471 

.000001 

.004412 

-.006571 

.01 

.020 

0.00 

-.004642 

.004 .6< 

.000251 

.005767 

-.005641 

-1.27 

.020 

10.00 

-.001419 

.001412 

-.009991 

.004446 

-.006444 

18.21 

.020 

40.00 

.001910 

-.0n4o07 

-.000104 

.001762 

-.007443 

-12.02 

-  .010 

90.00 

.001841 

-.004116 

.001980 

.001949 

-.004442 

12.34 

-  .040 

0.00 

-.004690 

.004788 

.000417 

.004791 

-.004701 

-2.61 

.040 

10.00 

-.002844 

.OOlllO 

-.008420 

.001794 

-.004520 

IS.  14 

.040 

so.oa 

.001898 

-.004792 

-.006542 

.001211 

-.OOS124 

-11.46 

.040 

90.00 

.002844 

-.001118 

.001847 

.002986 

-.001349 

14.96 

DisPLACSiENTS  (Output  OptioD  No.  6) 


DIST 

ANOLB 

U 

V 

0.000 

0.00 

.000794 

-.000080 

0.000 

30.00 

.000941 

.000342 

0.000 

60.00 

.000691 

.000435 

0.000 

90.00 

.000214 

.000294 

.020 

0.00 

.000645 

-.000093 

.020 

30.00 

.000794 

.000300 

.020 

60.00 

.000418 

.000471 

.020 

90.00 

.000050 

.000234 

.040 

0.00 

.000418 

-.000104 

.040 

10.00 

.000668 

.000268 

.040 

60.00 

.000178 

.000419 

.040 

90.00 

-.000045 

.000192 

DO  YOU  WISH  TO  CONTINUE? 
7NO 

STOP 


OPIMIII-tSO 

Figur*  10.  Locd*d  Hol«  Ca«« 
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IDO  YOU  I/AIIT  INSTRUCTIONS? 

TYKS 

SGLCCT  DCSIRrn  OUTPUT  phoh  ths  polumino  casks. 
I  CARPET  plot  data 
LAMINATE  PROPERTISS 
LAMINATE  STRESSES 


2 

I 

4 

5 

6 

7 

8 
9 

10 
73. S 

INPUT  MUMDER  OF  DIFFERENT  PLIES  TO  BE  INPUT  (S  MAX)  AMD 
NUMBER  OF  DIFFERENT  MATERIALS  (I  MAX) 


LAMINATE  STRAINS 

CIRCUMFERENTIAL  4  RADIAL  STRESSES/ STRAINS 

DISPLACEMF.NTS 

STRAINS  PER  PLY 

STRESSES  PER  PLY 

FAILURE  CRITERIA  PF.R  PLY 

AUTOMATIC  SEARCH  FOR  FAILURE 


THE  UNIDIRCCTIONAI,  MATERIAL  PROPERTIES  FOR  BACH  MATBRXAL 
IN  THB  POLLOWXNG  ORDERt  Elp  E2,  012*  P0X880NS  RATXO 

-710,E€, 10.E6,3.B5E6. ,3  * - 

XNPUT  TtIB  ANQUIAR  ORIENTATION  Of  EACH  PLY 
70. 

Itlpirr  THE  THICKNESS  OF  EACH  PLY 


LT  I 

Lti 


Isotropic  Profjerties 
(e1  =  E2,  G= 


2(1  +1^) 


INPUT  I  FAR  FIELD  ST?.ESSE8  PX.PY.PXY,  OFF  AXIS  ANOLB,  BBARIHO  STRESS 
AND  BOLT  LOAD  1  NO  ANGLE. 

710000. . 0..0..0..0..0. 

INPUT  BOLT  DIAMETER,  DEGREES  BETWEEN  OUTPUT,  LOW  RjKNOE,  HIGH  RANGE, 
STEP  INCREMENT  AND  NUMBER  OF  STEPS  DESIRED  (7  MAX) 

7.25.19.. 0..90..0..1 


nisT 

o.ooo 

0.000 

0-000 

0.000 

0.000 

0.000 

0.000 


LAMINATE  STRESSES  (Output  Option  No.  3) 


ANOLB 

0.00 

15.00 

30.00 

45.00 

60.00 

75.00 

90.00- 


X  STRESS  Y  STRESS 


-.00 

-490.49 

-.41 

5001.63 

15003.65 

35480.07 

39993.51 


-10000.00 
-6031.64 
-1.33 
5001 .63 
5U01.33 
1830.03 
-.00 


SHEAR 
STRESS 
0.00 
ir  .0.53 
.70 
-50)1.63 
-86.  3.36 
-68;  9 . 73 
-.00 


MAX. 

MIN. 

DIRECTION 

PRINCIPAL 

PRINCIPAL 

-.00 

-10000.00 

0.00 

-.00 

-7322.1) 

15.00 

-.00 

-1.0.! 

30.00 

10003.25 

-.00 

45.00 

30004.87 

.OC 

-30.00 

273IB.88 

.00 

-15.00 

29993.51 

-.00 

-.00 

CIRCUMFCHENTIAL  AND  RADIAL  STRESSES  ii  STRAINS  (OutpUl  Option  No.  5) 


DIST 

ANGLE 

THETA 

RADIAL 

SHEAR 

THETA 

RADIAL 

SHEAR 

STRESS 

STRESS 

STRESS 

STRAIN 

STRAIN 

STRAIN 

0.000 

0.00 

-10000.00 

-.00 

0.00 

-.001000 

.000300 

0.000000 

0.000 

15.00 

-7322.13 

-.00 

-.00 

-.000733 

.000220 

.000706 

0.000 

30.00 

-1-62 

-.00 

,00 

-.000000 

.000000 

.000000 

0.000 

45.00 

10003.35 

-.00 

.00 

.001000 

-.000299 

-.000000 

0.000 

60.00 

20004.87 

.00 

.00 

.002000 

-.000599 

.000562 

0.000 

75.00 

27318.88 

.00 

.00 

.002731 

-.000019 

.000767 

0.000 

90.00 

29993.51 

-.00 

-.00 

.002999 

-.0009(50 

.000000 

DO  YOU 

WISH  TO 

CONTINUE? 

7N0 


STOP 


QP1 14)1 16-221 

Figura  11.  Isotropic  Unloadsd  Hols 
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IDO  YOU  WWfT  IMSTROCTIOMS7  ' 

7YES 

SELECT  DESIRED  OUTPUT  FROM  THE  EOLLOMINO  CASES.  f 

1  CARPET  PLOT  DATA  ^ 

2  LA’IINATE  PROPERTIES  \l 

3  LAtllHATE  STRESSES  ^ 

4  LAMIHATE  STRAIUS  1  . 

5  CiaCUHFEREUTIAL  «■  RADIAL  STRESSES/ STRAINS  ^ 

6  DISPLACEMENTS  ^  \ 

7  STRAINS  PER  PLY  'I 

B  STRESSES  PER  PLY  X  90°  j 

9  FAILURE  CRITERIA  PER  PLY  \  1  .  j/ 

10  AUTOMATIC  SEARCH  FOR  FAILURE  ^ 

710  \ 

INPUT  NUMBER  OF  DIFFERENT  Pt.IF.S  TO  DE  INPUT  (8  MAX)  \ 

AND  NUMBER  OF  DIFFERENT  MATERIALS  {i  MAX) 

' 

iNPOT  THE  UNIDIRECTIONAL  MATERIAL  PROPERTIES  FOR  EACH 
MATERIAL  IN  THE  FOLIOUING  ORDERi  El,  F,2,  012,  POISSONS  RATIO 
718. OSES, 1 .9E6,  .aSEO,  .3 

INPUT  THE  UNIDIRECTIONAL  Al.lXXJABLES  FOR  EACH  MATERIAL 
IN  THE  FOLLOUINO  OROERl  Tl.  Cl,  T2 .  C2,  SHEAR 
72.3E5.3.2E5,2.n2E4,  23E4,1.73E4 

INPUT  THE  ANGULAR  ORIENTATION  OF  EAai  PLY 

70. . 45. ,-45, ,90. 

INPUT  THE  THICKNESS  OF  EAOl  PLY 

7.5. . 2. .2,  .1 

INPUTt  FAR  FIKLO  STRESSES  PX.PY.FXY,  OFF  AXIS  ANGLE,  BEARING  STRESS 
AND  dolt  LOADING  ANGLE. 

710000 . .  500 . . 2500 . . 30 . . 35000 . . 10 , - 

INPUT  UIDTH  (0.0  FOR  INFINI-^E  PLATE) 

70.0 

INPUT  BOLT  DIAMETER,  DEGREES  DETWEEN  OUTPUT,  LOW  RANGE,  HIGH  RANGE, 
STEP  INCREMENT  AND  NUMBER  OF  STEPS  DESIRED  (7  MAX) 

?.25,5..  0.,3S0., .02 , 2 

lUPUr  file  NUHDER  WHICH  CORRESrOlJDS  TO  THE  FAILURK  CRITERIA 
YOU  WISH  TO  USE 

1  M.\xiMua  strain 

2  MAXIMUM  STRESS 

3  TSAl-niLL 

4  MODIFIED  TSAI-WU 

5  HOFFMAN 


•Input  Stresses  are  Ratioed 
Uniformly  Until  First  Ply 
Failure  is  Predicted 


AUTOMATIC  SEARCH  FOR  FAILURE  (OutpUt  Option  No.  10) 
FAILURE  STRESSES 

PX  PY  PXY  P  -| 

17243.40  842.17  4310.85  60351.80  - 


DIST 

AtIGLE 

k*LY 

5*M  LURE 

FMLURE  RATIOS 

NUMBER 

2 

020 

335,00 

0.00 

.355 

•.210 

.46') 

.05  3 

020 

335.00 

45.00 

.999 

.011 

-.101 

.537 

020 

335.00 

-45.00 

.730 

•.400 

.650 

.620 

020 

1 

335.00 

i 

90.00  — 

1 — l.OOO 

1  .606 

-.075 

aitj 

—  Indicates  Ply  Failure 

* — Failure  Initiation  Location  Angle 

User  Specified  Distance  Away  from  Hole 
at  Which  Failure  is  Predicted 


-Depicts  Relative  Magnitude  of 
Stress  Components  to  Failure  Number 


Note;  The  "search''  is  only  done  within  the  user  specified  "range'  with  an  accuracy  i.t  1/2  the 
"degrees  between  output"  in  locating  failure  initiation. 


GPISal  16-319 


Figure  12.  General  Loading  Condition 
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PROGRAM  LISTING 


o  o 


c 

C  BOLTED  JOINT  STRESS  FIELD  MODEL  (BJSFM) 

C 

C  THIS  PROGRAM  COMPUTES  LAMINATE  STRESS  AND  STRAIN  DISTRIBUTIONS 

C  AROUND  LOADED  AND  UNLOADED  FASTENER  MOLES.  THE  PRINCIPLE  OF 

C  SUPERPOSITION  IS  USED  TO  OBTAIN  STRESS/ STRAIN  DISTRIBUTIONS 

C  FOR  A  GENERAL  APPLIED  LOADING.  OPTIONAL  FAILURE  ANALYSIS 

C  ROUTINES  ALLOVJS  LA-MINATE  STRENGTH  PREDICTIONS  USING  VARIOUS 

C  MATERIAL  FAILURE  CRITERIA  AND  HYPOTHESES.  PROGRAM  '/AS 

C  DEVELOPED  BY  J  .M.OGONOl/SKI  OF  MCDONNELL  AIRCRAFT  CO., 

C  ST. LOUIS,  MISSOURI, 

C 

c 

c 

PROGRAM  BJSFM( INPUT=100,OUTPUT=150,TAPE105»INPUT, TAPE5»INPUT, 
1TAPE10Q=OUT PUT. TAPES) 

C 

COMMON/  0;iE/El(3),E2(3),G12(3),V12(3) 

COMMON/ TWO/ lOUT (15), NUMPLY, NUMMAT, ANG ( 9 ) , PLYTHK { 3 ) , MATID( 9 ) 
COMflON/THREE/IANG,  II/IW ,  IHIGH,  STPI NK,  NUMSTP 
COMMON/ FOUR/ PX, PY.PXY, P, PW, ALPHA, BETA, DIA, CORRECT 
COMMON/ FIVE/ FXT ( 3 ) , FXC ( 3 ) , FYT( 3 ) , FYC ( 3 ) , FXY ( 3 ) , I FAIL 
COMMON/SIX/AI(3, 3) 

COM MON/ SEVEN/ S( 3, 3) 

COMMON/ EIGHT/ STRESS( 3 ,7,73), STRAIN( 3, 7 , 73 ) 
COMMON/NINE/STRl(B,7, 73) ,STR2(B, 7,73) , STRl 2 { 3 , 7 , 7 3 ) 

INTEGER  ANS,ANS2,YE3 
REAL  lANG,  ILOV/.IHIOH 
C 

DATA  YES/ ’YES' /.NO/ 'MO*/ 


10  CONTINUE 

OUTPUT(6) ' IDO  YOU  WANT  INSTRUCTIONS?' 

READ(5,20)ANP 
20  FORMAT (A3) 

C 

IF(  ANS.EQ.YES)  '/RITE(6,30) 

30  FORMAT (  *  SELECT  DESIRED  OUTPUT  FROM  THE  FOLLOl/ING  CASES.’/, 

1'  1  CARPET  PLOT  DATA’/’  2  LAMINATE  PROPERTIES’/’  3  ', 

2 'LAMINATE  STRESSES’/’  4  LAMINATE  STRAINS ’/ ’  5  CIRCUM*:”  , 

3’ERENTIAL  S<  RADIAL  STRESSES/ STRAI  NS’ /  ’  6  DISPLACEMENTS’/, 

4’  7  STRAINS  PER  PLY’/’  8  STRESSES  PER  PLY’/, 

5'  0  FAILURE  CRITERIA  PER  PLY’/’  10  AUTOMATIC  SEARCH  FOR  ’, 

6 ’FAILURE’ ) 

C 

DO  40  Ii=l ,  15 

IOUT(L)-0 

40  CONTir/UE 

REAn(5,  ) (lOUT(L) ,L=l , 10) 

C 

IF(  ANS.  EQ.  YES)  UP,ITE(o,  50) 

50  F0IV1AT(  ’  INPUT  NUMBER  OF  DIFFERENT  PI.IES  TO  BE  INPUT  (8  MAX)  ’  , 
1’  AND’/'  NUMBER  OF  DIFFERENT  MATERIALS  (3  MAX)') 

READ(5,  ) NUMPLY, NUMMAT 

C 

IF(  ANS.  EQ.  YF.S)  WRITE(  5,  GO) 

60  FORMAT( ’  INPUT  THE  UNIDIRECTIONAL  MATERIAL  PROPERTIES  FOR  ’, 
I’EACH  MATERIAL’/’  IN  THE  FOLLO'NING  ORDER;  El,  E2 ,  G12,  ’, 

2 'POISSONS  RATIO’) 

READ(5,  ) (E1(L) ,E2(L) , G1 ? ( L) , VI 2 ( D . L=1 , NUMMAT) 

C 
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C 

TFiANS.BQ.YES.AND.PUTOUT(IOUT,9) .GB.l. )  WRITE(6,70) 

70  PORMAT( '  INPUT  THE  UNIDIRECTIONIO^  ALLOWABLES  POR  BACH  ' « 

IMMATERIAL'/’  IN  THE  POLLOWIMG  ORDERt  T1 ,  Cl,  T2,  C2,  SHEAR') 
IF(PUTOUT(IOUT.9) .GE. 1. )  READ(5,  ) (FXT( J ) , rXC( J ) , PYT( J) , 

IFYC ( J ) , rXY ( J ) , J«1 , NUMMAT) 

0 

INPUT  ANGULAR  ORIENTATION.  THICKNESS  AND  MATERIAL 
IDENTIFICATION  PER  PLY 

IF^ANS.EQ.YES)  WRITE(6,80) 

80  FORMAT V '  INPUT  THE  ANGULAR  ORIENTATION  OF  EACH  PLY’ ) 

READ(5,  ) {ANG(L),L-1,NUMPLY) 

IF(ANS.EQ.YES)  WRITE(6,90) 

90  FORMAT ( '  INPUT  THE  THICKNESS  OF  EACH  PLY' ) 

READ(5,  ) (PLYTHK(J),J-1,NUMPLY) 

100  CONTINUE 

IF(NUMMAT.NE. 1)  GO  TO  120 
DO  110  L-1,NUMPLY 
MATID(L)-1 
110  CONTINUE 
GO  TO  140 
120  CONTINUE 

IF(ANS.EQ.YES)  WRITE{6,130) 

130  FORMATC  INPUT  THE  MATERIAL  NUMBER  FOR  EACH  PLY’) 

READ(5,  ) (MATID(L) ,L-1,NUMPLY) 

140  CONTINUE 

IF(PUTOUT(IOUT, 3) .EQ.0.0)  GO  TO  200 

IF'.  ANS.EQ.YES.  AND.PUTOUT(  lOUT,  3)  .GE.  1. )  WRITB(6, 150) 

150  FORMATC  INPUT:  FAR  FIELD  STRESSES  PX.PY.PXY,  OFF  AXIS  ANGLE,’, 

1’  BEARING  STRESS’/’  AND  BOLT  LOADING  ANGLE.’) 

IF(PUTOUT(IOUT.3)  .GE.l.)  READ{5,  )  PX,  PY,  PXY,  BETA,  P,  AI.PHA 

IF(ANS.EQ.YES.AND.P.NE.O.O)  WRITE(6, 160) 

160  FORMATC  INPUT  WIDTH  (0.0  FOR  INFINITE  PLATE)’) 

IF(PWrOUT(IOOT,  3)  .GE.  1.  .AtID.P.NE.O.O)  READ(5,  )  W 
C 

170  CONTINUE 

IF(ANS.EQ.YES. AND.PUTOUT( lOUT, 3) .GE.l.)  WRITE(6, 190) 

180  FORriAT(  ’  INPUT  BOLT  DIAMETER,  DEGREES  BETWEEN  OUTPUT,  LOW  ’  , 

1 'RANGE,  HIGH  RANGE,'/’  STEP  INCREMENT  AND  NUMBER  OF  STEPS’, 

2’  DESIRED  (7  MAX)’) 

C 

IF(PUTOUT( lOUT, 3) .GE. 1. )  READ(5,  )  DIA, lANG, ILOW, IHIGH, 
ISTPINK.NUMSTP 
C 

IF(NUMSTP.GT-- 7)  NUMSTP-7 
RANGE«(  IHIGH-ILaO/IANG 

IF( RANGE. GT. 72)  OUTPUT(6)’  DEGREES  BETWEEN  OUTPUT  TO  SMALL;  ’, 

1’  TRY  AGAIN’ 

IF(RANGE.GT. 72)  GO  TO  170 
C 

BL-P*DIA 
PW-0 . 0 

IF(P.NE.O.O.AND.W.NE.O.O)  PW-BL/ ( 2 . 0*W) 

WD-W/DIA 

IF(WD.LT.4.0.AND.VJ.NE.0.0)  OUTPUT(6)’  CAUTION:  WIDTH-TO-DIAMETER’ 
1'  RATIOS  LESS  THAN  4  GIVE  ERRONEOUS  RESULTS’ 
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IF(AIir;.  EQ.  VT:S.ANn.PUTOUT(IOUT,  9)  .GE.  1  .  )  WRITE(6, 190) 

190  FORMAT!'  INPUT  THE  NUMBER  l/llICU  CORRESPONDS  TO  T!iE 

1’ FAILURE  CRITERIA  '/*  YOU  WISH  TO  USE*/*  1  MAXIMUM  STRAIN*/, 
2*  2  MAXIMUM  STRESS*/*  3  TSAI-HILL*/*  4  MODIFIED  *, 
3*TSAI-V/U*/*  5  HOFFMAN*) 

IF{PUTOUT(IOUT,  9)  .GE.  1.  )  READ(5,  )IFML 

IF(PUTOUT!IOUT,  10)  .EQ.2.  )  NU>4STP*2 
200  CONTINUE 

CARPET  PLOT  ROUTINE 

IF(PUTOUT(IOUT, 1) .NE.2. )  GO  TO  240 
A!IG1=ANG(  1  ) 

ANG2“ANG( 2) 

ANG3»ANG(3) 

ANG4=ANG ( 4 ) 

FLYT1=*PLYTHK(  1 ) 

PLYT2=PLYTHK(2) 

PLYT3=PLYTMK(3) 

FLYT4=PLYTHK(4) 

DO  200  JKI=1. 11 
FORFIV=-. 1 
DO  270  IJK=1,11 
ANG( 1 )=ANG1 
ANG( 2)=ANG2 
ANG( 3 )=ANG3 
AI1G(4)=ANG4 
PLYTHKd  )  =  PLYT1 
PL\'TH:'(2)=PLYT2 
PLYT1IK(3)=PLYT3 
PLYT‘IK(4)=PLYT4 
CJUMPLY=4 

FLYTHK! 1 )»1 . 1-JKI*  .  1 
FORPTV=<FORFIV+  .  I 
CHECK=*FORFIV+PLYTHK(  1 ) 

1F(CI1ECK.GT.  1 . 0)  GO  TO  280 
PLYTUK! 2 )=FORFIV/2 . 

PLYTHK ( 3 ) »FORFI V/ 2 . 

PLYTHK(4)=*1 . 0-PLYT'lK(  1  )-FORFIV 
I0=PLYT'.1K(  1  )  *100.4 
I45=PLYTHK( 2) *200.4 
I90=PLYTHK(4)*100.4 
IF(IO.NE.O)  GO  TO  210 
NUMPLY-NUMPLY-1 
PLYTHKd  )=»PLYTHK(2) 

PLYTHK(2 )=PLYTHK(3 ) 

PLYTHK ( 3 ) =PLYTaK ( 4 ) 

ANG!1)=ANG{2) 

A:1G(2)=ANG!  3) 

ANG(3)»ANG(4) 

210  CONTINUE 

IF(I45.NE.O)  GO  TO  220 
N'JMPLY=NUMPLY-  2 

IF(NUMPLY.EO. 1 )  PLYTHK {1 ) =PLYTHK{ 4 ) 

IF(NUMPLY.EQ.  1 )  AHGd)=ANG{4) 

PLYTnK(2 )=PLYTHK(4) 

ANG(  2)=At.’G(4) 

220  CONTINUE 

IF (  1 90 .  EQ .  0  )  NUMPLY»NL'MPLY-1 

IF(PUTOUT(IOUT, 1) .EQ.2. )  WRITE(6,^30)  10,145,190 
230  FORt4AT(//25X, 'LAYUP;  *  ,  13  ,  *  /  *  ,  1 3  ,  *  /  *  ,  1 3  ) 
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BRANCHES  TO  OTHER  SUBROUTINES 

240  IF(PUTOUT(IOUT, 2) .EQ.O. )  GO  TO  260 
ALPH-ALPHA 
CALL  ABD(ALPH) 

CORRECT-1 . 0 
OUMMY-PUTOUT{ lOUT, 98) 

250  CONTINUE 

IP(PUTOUT(IOUT, 3).EQ.0.)  GO  TO  260 
CALL  LAMSTR 

1F(PUTOUT(IOUT,7).EQ.O.)  GO  TO  260 
CALL  PLYSTR(IFAIL) 

IF(PUTOUT(IOUT,9) .EQ.O. )  GO  TO  260 
CALL  FAILURE 

IF(PUTOUT{IOUT, 10) .EQ.2.)  DUMMY-PUTOUT( lOUT, 99) 
IF(CORRECT.LT. .999. OR. CORRECT. GT. 1.001)  GO  TO  250 


260  CONTINUE 

DUMMY-PUTOUT  (  I  OUT  ,.90 ) 

IF(PUTOUT(IOUT, 1) .NE.2.)  GO  TO  290 

270  CONTINUE 
280  CONTINUE 
290  CONTINUE 

OUTPUT(6)'  DO  YOU  WISH  TO  CONTINUE?’ 
READ(5,300)  ANS2 
300  F0RMAT(A3) 

IF (ANS2.NE.no)  go  to  10 


STOP 

END 

********************************************************* 

FUNCTION  PUTOUTdOUT,  IN) 

DIMENSION  IOUT(15) 

PUTOUT-0 . 

DO  10  J-1,15 

IF(IOUT(J) .GE. IN)  PUTOUT=l. 

10  CONTINUE 

DO  20  J-1,15 

lF{IOUT(J) .EQ.IN)  PUTOUT»2. 

20  CONTINUE 

IF(IN.EQ.98)  DATA=0.0 
IF(IN.EQ.99)  DATA-1.0 

IF ( DATA . EQ . 1 . 0 . AND . I N . LT . 1 0 )  PUTOUT-0 . 5 

RETURN 

END 
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SUBROUTINE  ABD  (ALPtlA) 

CALCULATES  LAMINATE  PROPERTIES  PROM  LAMINA  INPUT 
CONMON/ONE/El(3)  .E2(3)  .0120)  ,V12(i) 

COMMON/TWO/IOUT(15) ,NUMPLY,NUMMAT.ANG(8) ,PI  YTHK(8) ,MATID(8) 
COMtiOH/SIX/AI(3,3) 

COMMON/S EVEN/S( 3, 3) 


DIMENSION  V21(  3)  .DTVO)  .011(3)  ,022(3)  .012(3)  ,066(3)  ,U1( 3)  , 
1U2(3) ,U3(3) .U4(3) ,U5 ( 3 ) .0BAR(8, 3. 3 ) . ZZ( 16) . Z(  16)  ,0(3. 3) , 
2AA(3.3),A(3.3) 


REDUCED  STIFFNESSES  FOR  EACH  MATERIAL 

DO  10  M-1,NUMMAT 
V21(M)-E2(M)*V12(M)/E1(M) 

DIV(M)-1 .-V12(M)*V21(M) 

011(M)-E1(M) /DIV(M) 
022(M)-E2(M)/0IV(M) 
012(M)-V12(M)*E2(M)/DIV(M) 
066(M)-G12(M) 

10  CONTINUE 

INVARIANT  PROPERTIES 


DO  20  M»1,NUMMAT 

U1(M)»(3.*011(M)+3.*022(M)+2.*Q12(M)+4.*Q66(M) )/8. 
U2(M)-(01l(M)-Q22(M) )/2. 

U3(M)«(011(M)+022(M)-2.*012(M)-4.*066(M) )/8. 
U4(M)-(011(M)+Q22(M)+6.*Q12(M)-4.*066(M) )/8. 
U5(M)«(011(M)+022(M)-2.*012(M)+4,*066(M) )/8. 

20  CONTINUE 

DO  30  I»l,3 
DO  30  J»l, 3 
A(I,J)=0. 

AA(I, J)»0. 

30  CONTINUE 

TRANSFORMED  REDUCED  STIFFNESSES  PER  PLY 
THICK*0. 

DO  40  L=1,NUMPLY 

DEG-ANa(L) •3.1415926535/180.0 

M»MATID(L) 

QBAR(L, 1 , 1 )»U1 (M) +U2(M) *COS( 2 .*DEG )+U3 (M) *COS(4 . *DEG) 
QBAR(L,l,2)=U4(M)-U3(M)*COS(4.*DEG) 

QBAR(L,  2, 2)=»Ul(M)-U2(M)*COS(2.*DEG)+U3(M)*COS(4.*DEG) 
QBAR(L,  1,3)=«.5*02(M)*SIN(2.*DEG)+U3(M)*SIN(4.*DEG) 
QBAR(L,2, 3)*.5*U2(M)*SIN(2.*DEG)-U3(M)'*SIN(4.*DEG) 
QBAR(L, 3, 3)»U5(M)-U3(M)*COS(4.*DEG) 

QBAR(L,  2, 1)=*QBAR(L,  1,2) 

QBAR(L,3,1)»QBAR(L,1,3) 

QBAR(L, 3,2)=QBAR(L,2,3) 

THICK“PLyTHK(L)+THICK 
ZZ(L+1) -THICK 
40  CONTIIWE 

2(1)  — 1. •THICK/2.0 
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CALCULATE  A  I-IATRIX 

DO  70  1-1,3 
DO  60  J-I,3 
DO  50  L-1,NUMPLY 
Z(L+l)-Z(l)-»-ZZ{L+l) 

ZA^Z(L-t-l)-Z(L) 

A(I,  J)-A(I.  J)+QBAR(L,1,  .J)*ZA 
50  CONTINUE 

MATRIX  Q  AND  OQ  ARE  DUMMY  MATRICIES  USED  IN  CAI^CUTJVTIONS  INVOLVING 
THE  MANIPULATION  OF  OTHER  MATRICIES 

Q(I, J)-A(I, J)/THICK 
60  CONTINUE 
70  CONTINUE 

COMPUTE  A/THICK  INVERSE  MATRIX 
I STEP- 1 

CALL  INVERSE(Q,AI) 

LAMINATE  MID-PLAUE  PROPERTIES 

EXl-l.0/AI(l,l) 

EY1-1.0/AI(2,2) 

VXYl— EXl*AI(l,2) 

GXYl-1 ,0/AI{3, 3) 

IF{PUTOUT(IOUT,2) .EQ.2.)  WRITE(6,R0)  SXl, EYl.GXYl.VXYl 
80  FORMAT(/25X, 'LAMINATE  PROPERTIES* / *  EX  -  * , E9. 3, 2X, ' EY  - 
1E9.3,2X, 'GXY  -  ' , E9. 3, 2X, ' VXY  -  •,F5.3) 

CALCULATE  MATERIAL  PROPERTIES  FOR  OFF-AXIS  BOLT  LOAD 

TRANSFORMED  REDUCED  STIFFNESSES  PER  PLY 

THICK-0. 

ALPHA-ALPHA*  3 . 141 5926 53  5/ no .  0 
DO  90  L-1,NUHPLY 
DEG-kNG(L)*3. 1415926535/no. 0 
DEG-DEG- ALPHA 
M-MATID(L) 

QBAR(L, I, 1)-U1 (M)+U2(M)*C0S(2.*DEG)+U3(M)*C05(4.*Dna) 

QBAR(L, l,2)-U4(M)-U3(M)*COS(4.*DEG) 

QnAR(L, 2,2)-Ul(M)-U2(M)*COS(2.*DEG)+U3(M)*COS(4.*DCG) 

QBAR(L,  1,3)-.5*U2(M)*SIN(2.*DEG)+U3(M)*SIN(4.*DEG) 
QBAR(L,2,3)-.5*U2(M)*SIN(2.*DCG)-U3(M)*SIN(4.*DEG) 

QBAR( L,  3 »  3 )-U5 ( M) -U3 ( M) *C0S( 4 . ‘DEG ) 

QBAR(L, 2. 1)-QBAR(L, 1, 2) 

QBAR( L, 3 . 1 ) -QBAR( L, 1.3) 

QBAR(L,  3, 2)-QBAR{L, 2, 3) 

THICK-PLYTHK ( L ) +THICK 
ZZ(L+1)-THICK 
90  CONTINUE 

Z(l)— 1  .*THICK/2.0 
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CALCULATK  AA  IIATRIX 

no  120  1-1,3 
DO  110  J-1 , 3 
DO  100  L-1,MUMPLY 
'/(T.+  l  )-•».(  1  )  +  7.7.(T,4.1  ) 

ZA-Z(L+1 )-Z(L) 

AA{  I,.J)-AA(  I,  J)+QnAR(L,  I.  J)*ZA 

100  COOTI  IUK 

Q{I,  J)-AA(I,  J)/T11ICK 
110  COOTIWi: 

120  C0:1TI  JUK 

COriPlJTV:  AA/TMICK  IiiV3RSE  !1ATRIX 
1 STEP- 4 

CM^L  I^IVERSE(Q,S) 

OFF-AXIS  LA.MI'JATE  PHOPCRTIER 

EX2-1 .0/S( 1, 1) 

EY2-1 .0/S(2, 2) 

VXY2— EX2*S(  1,2) 

CXY2-1 .0/S( 3, 3) 


ui;tur!i 

END 

********••**«*•**•*****•«•***********• 

su'iROiJTi  JE  i*iv:;rs;:  (x.xi) 

CAECUUAT'CS  T'li:  IWERSE  OF  A  3X3  HArUIX 

ni  lF,NSIOU  X(  3,  3)  ,  'CI(3,3) 

COr  iMOO  I  STEP 


D::T»(X(  1  ,  1  )*X(  2, 2)*X(  3. 3)  )  +  {X(l ,  2)  *X(  2,  3)*X(  3, 1 )  )  + 

1  (X(1,3)*X(2.1)*X(3,2))-(X(1,3)*X(2. 2)*X(3,1) )- 

2  (X(l,  1)*X(2,  3)*X(3,2)  )-(X{l,2)*>:(?,  1)*X(3,  3)  ) 

IF(DF,T.nQ.0.0)  GO  TO  10 

XI  (  1 , 1 ) - ( X ( 2 . 2 ) *X ( 3 , 3 ) -X ( ? , 3 ) *X ( 3 . 2 ) ) /OCT 
XHl,2)-(X(2,  3)«X(3,l)-XC’,l)*x(3,  3)  )/ni/r 
XHI,3)=>(X(2,1)«X(3,?)-X(2,2)*X(3,  1)  )/:)i:t 
XI  ( 2 , 2 )= ( N ( 1 , 1 ) *X ( 3 , 3 ) -X( 1 , 3 ) *X ( 3 , 1 ) ) / DFT 
XI(2,3)  =  (X(1,2)*X(3,1)-X{1,1)*X(3,2)  )/D1:T 
XI(3,3)-(X(l,l)*X(2,2)-X(l, 2)*X(2,1))/DET 
XI(2. 1)-(X(3, 2)*X(1, 3)-Xvl,2)*X(3, 3) )/DCT 
XI(  3,l)-(X(l,  2)*X(2, 3)-X{2,2)*X(l,  3)  )/D1:T 
:CI(3, 2)-(X(2, 1  )*X(l,  3)-X(l,  i)*A(2, 3)  )/DET 
GO  TO  30 

n  PRI''IT  20,I:'TEP 

20  rOR  SAT(  ■  SIJ  IROUTI  JE  lr5V:;,(SE  CAI,CWLATES  A  SINGULAR  MATRIX 
1  •  AT  STFJ*  '13,' 

30  CO’JTIMUE 
RETU'FJ 
E'lD 
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SUBROUTINE  LAMSTR 

CALCULATES  THE  TAMINATE  STRESSES  AND  STRAINS  DUE  TO  A 
OBNER^Ui  INPLANE  LOADING  WITH  A  DOLT  LOAD 

COMMO^  /TWO/ IOUT( IS). NUMPLT , NUMMAT. ANO( 9 ) , PLYTHK ( 9 ) , MATID(  9 ) 
OOMNON/THRBB/ lANG , ILOW . IHIGH.  STPINK, NUMSTP 
OONNDN^POUR/PX, PY. PXY, P, PW. ALPHA, BETA, DIA, CORRECT 
COMMON/SIX/AI(3,3) 

COMMON/ SEVEN/ S ( 3 .  J ) 

COMMON/ BIGHT/ STRESS^  3,7,73), STRAIN( 3,7,73) 

REAL  lA'IG.ILa),  IHIGH 

DIMENSION  STR(3,7,73),U(7.73),V(7,73) .UX( 7,73), VY( 7,73 ^ 

PX-CORRECT*PX 
py-correi:t*py 
PXY-CORRECT*  PXY 
P-C<»RECT*P 
PW^CORRECi'*PW 


PI-3. 14159;;  6;>35 
NUMPT-(  (IH3:GH-IL0W)/IANG)  +  1 
DO  10  J-1,NUMSTP 
DO  10  K-l,NliMPT 
U(J.K)-0.0 
V(J,K)-0.0 
DO  10  1-1,3 
STRESSd,  J,K)-0.0 
STRAINd,  J,K)-0.0 
10  CONTINUE 

CALCULATE  UNLOADED  HOLE  STRESSES 

li^(PX.EQ.O.O)  GO  TO  20 
BBTAO-BETA 

CALL  UNLODED(PX,DIA,AI, BETAO, STRESS, U,V) 
20  CONTINUE 

IF( Py.EQ.0.0)  GO  TO  40 
BETA90-BETA+90 . 0 

CALL  UNLODED ( P Y , D I A . AI , BETA90 , STR , UX , VY ) 
DO  30  J-1,NUMSTP 
DO  30  K-1,NUM1»T 
U(  J,K)-U(J,K)-tUX(J,K) 
V(J,K)-V(J,K)4VY(J,K) 

DO  30  1-1,3 

STRESS (I , J , K ) -STRESS (I , J , K )+STR( I , J , K ) 

30  CONTINUE 

40  CONTINUE 

IF("XY.EQ.0.0)  GO  TO  70 
BETA45-BETA+45 . 0 

CALL  UNLODED  (  PXl',  DIA,  AI ,  BETA45 ,  STR,  UX ,  VY ) 
DO  50  J~1,NUMSTF 
DO  50  K-1,NU?1PT 
U(J,K)-U(J,K)+UXtJ,K) 

V(J,K)-V(J,K)+VYt J,K) 

DO  50  1-1,3 

STRESS (I , J , K ) -STRESS (I , J . K ) +STR  d  .  J . K ) 

50  CONTINUE 
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DET^4S■r^^;TA-45.0 

pxy’:=«-p>cy 

CAI.f.  !J  JUOUEO(PXYN,DIA,AI,  QETA45,STR,tJX,VY) 
DO  GO  J-1,MUMSTP 
DO  CO  'C-l,fJU'lPT 
U(J.K)>iU(J,K)  +  UX(J,K) 

V(J,K)-V(J,  IC)+VY(J.K) 

DO  CO  I-l, 3 

STRORSd,  J,  !0=»*;'"KEr.S{l,  J,K)4-STIl(I,  J,K) 

60  COOTINUL 

c 

70  COOTItJUC 

c;UjCulate  LOAonr  hole  stresses 

IF( P. EQ.O.O)  OO  TO  100 
<U,  ’)1A0«ALPHA 

PBtitp 

CALL  L0ADEI)(  PB,  DIA,  0,  ALPHAO,  STH,  UX,  VY) 

DO  no  J»1,WUMSTP 
DO  ao  K«l,NUrlPT 
U(J,K)-U( J, K)+UX(J,K) 

V(J,K)«V(J,K)+VY(J,K) 

DO  no  I-1,J 

STRESS { I , J . K ) « STRESS ( I ,  J , K ) +3TR { I ,  J ,  K ) 

00  CONTI  flUE 

ALPHAO- ALP!  I A 

CALL  UNLODnD(P\J,OIA,AI,  ALPHAO,  STR,UX,VY) 

DO  90  J-1,HUMSTP 
DO  90  K-l,OUMPT 
U(J,K)-U( J,K)+UX(J,K) 

V(  J,K)-V(J,K)+VY( J,K) 

DO  90  1-1,3 

STRESad,  J,  K)-STUE.OS(I,  J,K)+GTU(I,  J,K) 

90  CONTIIIUE 
100  COtITIHUE 

lElPUTOUTdOUT,  3)  .EQ.  2.  )  WRITE(G,  110) 

110  EOPx:lAT(///2nX,  '  I,AriINATE  STRESSES'//’  OIST 
1'  Y  STRESS  SHEAR  flAX.  MIN. 

239X,' STRESS  PRINCIPAL  PRINCIPAL’) 


ANGLE  X  STRESS', 
DIRECTION' /, 
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CAU:UIATE  PRINCIPAL  STRESSES 

IF(PUTOUT(IOUT,3).NE.2.)  GO  TO  140 
DO  130  JJ-l.NUMSTP 
DO  130  NN-1,NUMPT 

PRINA-(STRESS( l,aj,NN)-STRESS(2. JJ.NN) ) * ( STRESS( 1 , JJ,  NN) - 
1  STRESS(2. JJ.NH) )/4. 

?R INA-SQRT(  PRTNA-»-STRESS  (  3 ,  J J  ,  NN)  ‘STRESSC  3  ,  J J ,  NN )  ) 

PRINl-  ( STRESS  ( 1 ,  a  J ,  !1N  )  +r.  TRESS  (  2  ,  J  J ,  NN  )  )  /  2  .  ♦PR  INA 
PRIN2-  ( STRESS{  1 ,  J  J ,  NN)  +S  i'RESS(  2 ,  J  J ,  NN )  )  /  2  .  -PRINA 
TSTS-STRESS(  1,  JJ.NN)-STRE::3(2,JJ.NN) 

DIRCT-0. 

IF(  TSTS.NE.  0  . )  OIRCT-.  5*A1  AtJ(  2 .  *STRESS(  3  ,  JJ,  NN)  /TSTS) 
DIRCT-lflO.*DIRCT/3. 1415926535 

IF(PUTOUT(  lOUT,  3)  .NE.  2.)  GC>  TO  140 
ANGLE- ( NN- 1 ) * I ANG+ I LOW 
DIS1'-(JJ-1  )*STPtNK 

WRITE(6,  120)  DIST,  ANGLE,  STRF.SSC  1 ,  JJ,NN)  ,  STRESS (  2  ,  JJ ,  NN)  , 
lSTREt;S(  3,  JJ,NN)  ,  PRINl ,  PKIN2 ,  DIRCT 
120  FORMAT (F6.3,Fa.2,5F11.2,FQ.2) 

130  CONTINUE 
140  CONTINUE 

IF(PUT0UT( lOUT, 4) .E0.2. )  WRITE(6, 150) 

150  FORMAT( ///28X, *  LAMINATE  STRAINS*//'  DIST  ANGLE  X  STRAIN*, 
1*  Y  STRAIN  SHEAR.  M-\X.  MIN.  DIRECTION*/, 

239X,' STRAIN  PRINCIPAL  PRINCIPAL*) 

CALCULATE  LAMINATE  STRAINS 

DO  160  JJ-1,NUMSTP 
DO  160  NN-1,NUMPT 
DO  160  KK-1,3 
DO  160  MM-1 , 3 

STRAIN ( KK , J J , NN ) -AI ( KK, MM) ‘STRESS! MM, J J , NN) +STRAIN ( KK , J  J , NN ) 
160  CONTINUE 

CALCULATE  PRINCIPAL  STRAINS 

IF(PUTOUT(IOUT,4).NE.2.)  GO  TO  1)0 
DO  189  JJ-1,NUMSTP 
DO  ISO  NN-1,NUMPT 

PRINA- ( strain! 1 , JJ , NN) -STRAIN! 2 , JJ . NN ) ) * ! STRAIN ! 1 , J  J , NN ) - 
1  STRAIN!2,  JJ,NII)  )/4. 

PRI NA-SQRT  !  PR  I NA  I-  STRAI  N!3,JJ,NH)*.25*STRAIN13,JJ,Ntl)) 
PRIN1-!STRAIN!1, JJ,NN)+STRAIN!2, JJ,NN) )/2.+PRINA 
PRIN2-! STRAIN! 1 , JJ,NN)+STRAIN!2, JJ,NN) )/ 2. -PRINA 
TSTS-STRAI N !  1 ,  J  J ,  NN )  -SiTRAI  N  i  2 ,  J  J ,  NN  ) 

DIRCT-0 . 

IP ! TSTS . NE . 0 . )  DIRCT- . 5* ATAN! 2 . ‘STRAIN! 3 , J J , NN ) /TSTS ) 
DIRCT-180.‘DIRCT/3 . 1415926535 
D1ST-! JJ-1 ) ‘STPINK 

ANGLE-!  NN-1  )*IANG+ILO^■/ 

WRITE! 6 , 170 )  DIST,  -'VNGLE,  STRAIN!  I ,  JJ,  NN)  ,  STRAIN!  2  ,  JJ,  NN)  , 
1STRAIN!3. JJ.NN) , PRINl , PRIN2 , DIRCT 
170  FORMAT!F6.3,F8.2,5Fll.6,F8.2) 

180  CONTINUE 
190  CONTINUE 
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CALCUIATE  CIRCUMFERENTIAL  AMD  RADIAL  STRESSES  {■  STRAINS 


IF( PUTOUT( lOUT, 5 ) .EQ. 2. )  WR1TE(6. 200) 

200  FORMAT!// A  5X, 'CIRCUMFERENTIAL  AND  RADIAL  STRESSES  S,  STRAINS', 

1//'  DIST  -ANGLE  THETA  RADIAL  SHEAR  THETA', 

2'  RADIAL  SHEAR' /21X, 'STRESS  STRESS 

3'STRESS  STRAIN  STRAIN  STRAIN') 

IF!PUTOUT(IOUT,5).NE.2.)  GO  TO  230 
DO  220  J-1,NUMSTP 
DO  220  N-1,NUMPT 

ENERGY- . 5* ( STRESS ( 1 , J , N ) *STRAIN( 1 , J , N ) +STRESS( 2 , J ,  N ) * 

1  STRAIN(2, J,N)+STRESS(3, J,N)*STRAIN(3, J,N) ) 

ANGLE- (  N- 1  )  •  I ANG+ 1 LOV; 

D-ANGLE*PI/ia0  0 
DIST-( J-1)*STPINK 

RADSTS-STRESS ( I , J , N ) *COS ( D ) *COS( D ) +STRESS (2,J,N)*SIN(D)* 

I  SIN{D)+2.*STRESS{3, J,N)*SIN(D)*COS(D) 

CIRSTS-STREGSCI, J,N) *SIN!D)-SIN(D)+STRESS(2, J,N)*COS(D)* 

1  COS( D)-2. ‘STRESS! 3, J,N)*SIN(D)*COS(D) 

SHRSTS— 1 . *  STRESS! l.J,N)*SIN!D) ‘COS ! D ) +STRESS !2,J,N)*SIN(D)* 

1  COS!D)+STRESS!3,J,N)‘!COS!D)*COS!D)-SlN!D)‘SINfD) ) 

radstn-strain!i, j,n)‘cos!d)‘cos!d)+strain!2, j,N)*siN!n)‘ 

I  sin!d)+strain!3, j,n)‘sin!d)‘cos!d) 

CIRSTN-STRAIN!l;J,N)‘SIN!D)‘SIN!D)+STRAIN!2, J,U)*C0S!D)‘ 

1  C0S!D)-STRAIN!3, J,N)‘SIN!D)*COS!D) 

SHRSTN—1. ‘STRAIN!  1,J,N)‘SIN!D) ‘COS!  D)+STRAIN!  2,  J,N)‘SIN!D)* 

1  cos!n)+qTRAiN! 3, j,n)‘!cos!d)‘cos!d)-sin!d)*sin!d) ) 

\«RITE!6 , 210  )  DIST,  ANGLE ,  Cl RSTS ,  RADSTS,  SHRSTS,  CIRSTN,  RADSTN, 

1  SHRSTN 

210  FORMAT!F6.3,F10.2,F12.2,2F11.2,3F9,&) 

220  CONTINUE 
230  CONTINUE 

OUTPUT  DISPLACEMENTS 

U3-3.0‘DIA 

DISP-DIA/2 .0+NUMSTP‘STPINK 

IF ! PUTOUT ! I OUT ,  6 ) . EQ . 2 . 0 . AND. P.NE.O.O.AND.DISP.GT.D3)  OUTPUT ! 6 ) 

1'  CAUTION:  DISPLACEMENTS  AT  POINTS  GREATER  THAN  3D  AWAY  ', 

2'  FROM  THE  HOLE  MAY  OE  IN  ERROR' 

IF!PUTOUT!IOUT,6) ,EQ.2. )  WRITE!6,240) 

240  FORMAT! ///20X,' DISPLACEMENTS'//'  DIST  ANGLE' lOX, ‘ U ' lOX, 'V ) 

IF!PUTOUT!lOUT,6) .NE.2.)  GO  TO  270 
DO  260  J-1,NUMSTP 
DO  260  K»1,NUMPT 
ANGLE- !K-1 )‘IANG+ILOW 
DIST-! J-1 )‘STPINK 

WRITE! 6, 250)  DIST, ANGLE,U ! J, K ) . V! J , K ) 

250  FORMAT !F7.3,F10.2,F13.6,F12.6) 

260  CONTINUE 
270  CONTINUE 
RETURN 
END 
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c 

SUDROUTiNE  UNLODEDIP.DIA.AI. BETA, STRESS, U, V) 

CALCULATE  STRESS  DISTRIBUTION  AROUND  AN  UNLOADED  HOLE 

COMMON/THREE/ I ANG,  ILOVJ ,  IHIGH,  STPINK,  NUHSTP 
REAL  lANG, ILOW, IHIGH 

DIMENSION  STRESS(3,7,73),U(7,73),V(7,73),AI{3.3) 

DIMENSION  WORK(5) , COEF( 5 ) , RTR(4 ) , RTI ( 4 ) 

COMPLEX  R1 , R2 , COMPLX , XI 1 , XI 2 , COMl , COM2 , DENI , DEN2 , PHI 1 ,  PHI 2 
COMPLEX  Z,Zi,Z2,Pl,P2,Ql,Q2 

CALCULATE  COMPLEX  PARAMETERS 

INITIALIZE  COMPLEX  NUMBER:  S0RT(-1.) 

COMPLX* ( 0. , 1 . ) 

NUMCO-4 

COEF( 1)-AI(2, 2)*1000000 
COEF(  2)— 2.*AI(2, 3)*1000000 
C0EF(3)*(2.*AI( 1, 2)+AI(3, 3) )*I000000 
COEF(4)— 2.*  AI(  1,3)  *1000000 
COEF(  5)=»AI(  1,  1)  *1000000 
CALL  ROOTS (COEF,\roRK,NUMCO, RTR, RTI, IE) 

Rl*RTR(l )+COMPLX*RTl(l) 

IF(RTI(2) .GT.0.0)  F 1*RTR( 2 )+COMPLX*RTI { 2 ) 

R2=RTR( 3 )+COMPLX*RTI ( 3 ) 

IF(RTI(4) .GT.0.0)  R2-RTR(4)+COMPLX*RTI(4) 

Pl-AI (1,1) *R1 *R1+AI (1 , 2 ) -AI ( 1 , 3) *R1 
P2-AI (1,1) *R2*R2+AI ( 1 , 2 )-AI ( 1 , 3 ) *R2 
Q1-AI(1,2)*R1+AI(2, 2)/Rl-AI(2,3) 
a2=«AI(l,2)*R2+AI(2,2)/F;2-AI(2,3) 


PI*3. 1415926535 
BETA-BETA*  PI / 1 BO • 0 

NUUPT-(  (1HIGH-IL0V0/IANG)  +  1 

DC  20  JJ-l.NUMSTP 
DO  10  CItJ=l,NUMPT 

U(JJ,NN)-0,0 

V(JJ,NN)=0.0 

NNN-NM-1 

JJJ*JJ-1 

THETA'- :  UI»IANG+ILO\/)  *PI/130.0 
RADIUS=JJJ*STPINKt-DIA/2.0 

CALCULATE  X  fi.  Y  COORDINATES  OF  POINTS  AROUND  UNLOADED  HOLE 

X=RADIUS*COS( THETA) 

Y-RADIUS* SIN (THETA) 

CALCULATE  LOCATION  PARAMETERS  FOR  UNLOADED  HOLE  EQUATIONS 

Z1=X+R1*Y 

Z2=X-t-R2*Y 

Z=X+COMPLX*Y 


32 


nn  non  n  o  nnn  n  o  nno 


C 


MAPPING  FUNCTION 


XIl-CSQKT(Zl*Zl-DIA*DIA/4.-Rl*Rl*DIA*DIA/4. ) 
XI2-CSQaT(Z2*Z2-DIA*DIA/4.-R2*P2*DIA*DIA/4. ) 

CHOOSE  CORRECT  SIGN  OF  CSQRT 

XIl-Zl/XII 
XI2-Z2/XI2 

IF(REAL(XI1) .LT. -.00001)  XII— l.*XIl 
IF(REAL(XI2) .LT.-. 00001)  XI2— I.^XIZ 

XIl-1 .-XII 
XI2-1 .-XI2 

CALCULATE  PHI  PRIME 

CO^•ll-R2*SIN(  2  .  *BETA)  +  2  .  *COS(  BETA)  *C0S(  BETA)+CO^iPLX•  (  2  .  •R2« 
1  SIN(BETA)*SIN(BETA)+SIN(2.*BETA) ) 

C0M2-R1«SIN(  2.*DETA)  ♦■2  .*COS(QETA)  *COS(  BETA)+COMPLX*  (  2  .  *R1* 
1  3IN(UETA)*SIN(BETA)+SIN{2.*BETA) ) 

UEN1=*2  .  *0IA*  (  R1-R2  )  •(  1  .+COMPLX*Rl  ) 

DEN2=2 . *DIA* { R1-R2 ) * ( I . +COMPLX*R2 ) 

PHIl— C0MPLX*P*DIA*C0M1*XI1/(2.*DEN1 ) 
PHI2=COMPLX*P*DIA*COM2*XI2/(2.*DEN2) 

CALCULATE  STRESSES  AROUND  HOLE 

STRESSd,  JJ,NN)=«P*COS(DETA)*COS(BETA)+2.*REAL(R1*R1*PHII  + 

1  R2*R2*PHI2) 

STRESS(2.  JJ,MN)»P*SIN(UETA)*SIH(BETA)  +  2.*REAL(PHI1-»'PHI2) 
STRESSO,  JJ,NN)=<P*SIN(BRTA)*CCS(BETA)-2.*REALIR1*PHI1  + 

1  R2*PHI2) 

CALCULATE  DISPLACEMENTS 

XI1=1 .-XII 
XI2=^1  .-XI2 

XI1=Z1/XI1 
XI2=Z2/XI2 

DF,N1  =  16.*  (R1-R2)*(Z1+XI1) 

DEN2=16.* {R1-R2)*(Z2+XI2) 

PHI1=-P*DIA*DIA* (COMPLX+Rl )*COMl/DENl 
PHI2  =  P*DIA*DIA*  (C0MPT,X+R2)  *COM2/DFN2 

U(  JJ,NN)=2  .*REAL(  P1*PHI1+P2*PI1I2  ) 

V(  JJ,NN)=2  .  *REAL(''J1*PI1I1+Q2*PHI2) 

C 

10  CONTINUE 
20  CONTINUE 
C 
C 

RETURN 

END 
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SUBROUTINE  LOADED(P,DIA.,S, ALPHA, STRESS, U,V)  J 

CALCULATES  STRESS  DISTRIBUTION  AROUND  A  LOADED  HOLE  j 

ASSUMING  A  COSI'IE  BOLT  LOAD  DISTRIBUTION  ^ 

COMMON/ T\;0/IOUT(  15)  , NUMPLY, NUMMAT,  ANGO)  ,  PLYTHK(a)  ,MATIO(S)  1 

COMMON/TtlREE/IANG,ILO\^,IHIGH,  STPINK, NUMSTP  1 

H 

i 

REAL  lANG.ILOW.IHIGH  5 

COMPLEX  R1,R2,COMPLX,Z,Z1,Z2,CPOS(50),CNEG{50),CZERO,CM,  I 

lAKl ,  AK2 ,  XI 1 ,  XI2  ,  PHIl ,  P11I2 ,  GOMl ,  COM2 ,  XXI 1 ,  XXI 2  | 

COMPLEX  CHECK 1, CHECK 2, P1,P2,Q1,Q2 


COMPLEX  Al(50) ,A2(50) 

DIMENSION  AMATHX (4,4), DMATRX( 4 , 4 ) , STRESS{ 3 ,7,73) 
DIMENSION  U(7,73),V(7,73),S{3,3) 

DIMENSION  VraRK(5) , COEP{ 5 ) , RTR( 4 ) , RTI ( 4 ) 

INITIALIZE  COMPLEX  NUMBER:  SQRT(-1.) 

COMPLX= ( 0 . , 1 . ) 

CALCULATE  COMPLEX  PARAMETERS 
NUMCO=4 

COEF(  1)=S(2, 2)'^1000000 

COEF(2)=-2.*S(2, 3)*1000000 

COEF(3)  =  {2  .*S(l,  2)+S(3, 3)  )-'1000000 

COEF(4)=-2 .*S(1, 3)*1COOOOO 

COEFC 5 )=S( 1,1) *1000000 

CALL  ROOTS(COEF,vroRK,NUHCO, RTR, RTI, IE) 

R1=RTR( 1 )+COMPLX*RTI ( 1 ) 

IF(RTI(2) .GT.0.0)  Rl=RTR{2)+COMPLX*RTI(2) 

R2=RrR( 3 ) +COMPLX*RTI ( 3 ) 

IF(RTI(4) .GT.0.0)  R2=RTR(4)+COMPLX*RTI{4) 

P1=S(1, 1)*R1*R1+S(1,2)-S(1,3)*R1 
P2=»S(1,  1)*R2*R2+S(l,2)-S(1,3)*R2 
Q1=S( 1, 2)*R1+S(2, 2)/Rl-S(2, 3) 

Q2=S(lr  2)*R2  +  3(2,2)/R2-3(2, 3) 


PI=3. 1415926535 
THICK=0. 

DO  10  H=l, NUMPLY 
THICK=THICK+PLyTHK (H ) 
10  CONTINUE 
P*4.0*P/PI 


A  COSlrJE  LOAD  DISTRIBUTION  OVER  HALF  OF  HOLE  AT  AN  ANGLE 
ALPHA  TO  X  AXIS 

CALCULATE  COMPLEX  CONSTANTS 

PI2-PI/2.0 

M«-l 

20  CONTINUE 
M*M+1 

IF(M.EQ.l)  GO  TO  40 
30  CONTINUE 

C1-SIN( (M-1)*PI2)/(2*{M-1)) 

C2=*SIN(  (M+1)*PI2)/(2*(M+1)  ) 

C3»SIN{ (M-1 )*(-PI2) )/(2*(M-l) ) 

C4=SIH( (M+1)*(-PI2))/(2*(M+1)) 

C5=COS( (M-1 )*PI2)/ (2*(M-1) ) 

C6=*COS(  (M+1)*PI2)/(2*{M+1)  ) 

C7=»COS(  (M-1)*(-PI2)  )/(2*(M-l)  ) 

CQ=COS( (M+1)*(-PI2) )/(2*{M+l) ) 

CM=P*( (Cl+C2-C3-C4)-COMPLX*(-C5-C6+C7+C8) )/ (2.0*PI) 
IF'M.EQ.O)  CZERO=CM 
iF(M.GT.l)  CPOS(M)=CM 
IF(M.LT.-l)  MN=-1*M 
IF(M.LT.-l)  CMEG(MH)=CH 
IP{M.LE.O)  GO  TO  50 
M— l*M 
GO  TO  30 
40  CONTINUE 
C1=PI2 

C2=SIN(2.* (PI2) )/4. 

C3=SIN(2.*(-PI2) )/4. 

C4=SIN(PI2)*SIN(PI2)/2. 

C5=SIM(-Pl2)*SIU(-PI2)/2. 

CM=P*( (Cl+C2-C3)-M*COMPLX*(C4-C5) )/ (2.*PI) 

IF(M.EQ.l)  CPOS(l)=»CM 
IF(M.EQ.-l)  CNEG(1)=CM 
IF(M.EQ.-1 )  GO  TO  50 
M=-1*M 
GO  TO  40 
50  CONTINUE 
M=IABS(M) 

IF(M.LT.49)  GO  TO  20 

TRANSFORM  COMPLEX  PARAMETERS  INTO  REAL  AND  IMAGINARY  PARTS 


S1=REAL(R1 ) 
S2=RFAL(R2  ) 
T1=AIMAG( R1 ) 
T2=AIMAG(R2 ) 


C  EQUATING  COEFt’ICI EMTS  AND  SOLVING  FOR  CONSTANTS 

C 

DO  80  H-l,45 
MN-M-1 

IF(MN-NE.O)  GO  TO  60 
BMATRXd  )«REAL(-CZERO*DIA/2.  ) 

BHATRX(2)-AIMAO(-CZERO*OIA/2. ) 

GO  TO  70 
60  CONTINUE 

BMATRX(l)-REAL(-CPOS(MN)*DIA/(2.*{MM+l ) ) ) 
BtlATRX(2)-AIMAG(-CPOS(MN)*DIA/(2.*(MN+l)  )  ) 

70  CONTINUE 
MN-M+1 
MNEG— 1*^1 

Bt4ATRX(3)-REA;j(-CNEG(MN)*DIA/(2.*(MNEG+l)  )  ) 
BMATIOC(4)-AIKAG(-CNEG(MN)*DIA/{2.*(MNEG-H  ) ) ) 

AMATRXd,  1)-T1+1  . 

AMATRXd,  2)-Sl 
AMATRXd,3)»T2+l. 

AMATRXd,4)-S2 
At1ATKX(2, 1  )«S1 
A:1ATRX(2,2)»-T1-1. 

AriATRX(2, 3)-S2 
A:<ATRX(2,4)=-T2-1  . 

AMATRX(3. 1)-1 .-Tl 

AMATR:-(3, 2)  — SI 

AMATKX(3, 3)-l .-T2 

A'1ATRX(3,4J=— S2 

AJ1ATRX(4,  1  )-Sl 

AMATRX(4, 2)-»l  .-TI 

AflATRX(4, 3)-S2 

AilATRX(4,4)-l  .-T2 

CALL  SIMULT(AHATRX,B;1ATHX,4,J) 

IF(J.EQ.l)OUTPUT(6) '  SI  lULT  CALCULATES  A  SINGULAR  SET  OF  EQS 
A1  (M) =BMATRX( I )+COMPLX*nMATRX(  2 ) 
a2  ( M )  -UMATRX  (  3  )  '■COf^PLX*nMATRX(  4 ) 

80  CONTINUE 
C 

PX=2.*PI*AIMAG(COMPLX*DIA*CNEGd)/2.  ) 
PY-2.*PI*REAL(COMPLX*OIA*CNEGd  )/2.  ) 

C 

AIlATRXd,  l)»Tl 
AMATRXd,  2)-Sl 
AflATRXd,  3)»T2 
AMATRXd,4)-S2 
Ai'lATRX{2, 1  )-0 .0 

a;iatux(2, 2)=l  .0 

AMATRX(2,3)»0.0 
AflATRX(2,4)  =  l  .0 
AMATRX(3,1)=*2.»'51*T1 
AriATRX(3,2  )=Sl*r.l-Tl*Tl 
AMATRX(3, 3)=-2  .*S2*T2 
A:IATRX(3,4)«S2*52-T2«T2 
AMATRX(4, 1 )»- T1/(51*S1+T1*T1) 

A1ATRX(4,2)  =  S1/(51*S1+T1*TI  ) 

AMATRX(4,  3)=*-T2/(52*S2+T2*r?) 

AilATRX(4, 4)=S2/  (S2*S2+T2*T2) 

QMATRXd  )  =  PX/  (  4.*PI) 

3MATRX ( 2 ) =-Py/ ( 4 . *  PI ) 

BHATRX(3)-(S(1,2)*PY+S(l,  3)*PX)/(4.*PI*Sd,  1)  ) 
BHATRX14)--(S(1, 2)«PX+S(2, 3J*PY)/t4.*PI*S(2, 2) ) 

CALL  SINULT(AMATRX,DtlATRX,4,J) 

IF( J.EQ. 1 )OUTPUT(6) ■  SIMULT  CALCULATES  A  SINGULAR  SET  OF  EQS 
C 

AK  l«rj:iATRX  ( 1 )  +COMPLX*DMATRX  (  2 ) 

AK2=Bf UVTRX  (  3  )  +COMPLX*BMATRX  (  4 ) 
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C 


N UMPT-  (  ( I H IGH- 1 LOVJ )  / 1 ANG )  + 1 
Al..Pn A— ALPHA*  r  I / 1  BO .  0 
Al.PH— ALPHA 
DO  150  JJ«1,IIUHSTP 
DO  140  NN=«l,tIUMPT 

U(JJ.HN)»0.0 

V(JJ,NW)=0.0 

N:iU=NtI-l 

JJJ=JJ-1 

THET  A=  (  NUN*  1 ANG+ 1 LOVJ  )*PI/180.0 
RADIUS=J JJ*STPI WK+DIA/ 2 . 0 

CALCULATE  X  AND  Y  COORDINATES  OF  POINTS  AROUND  LOADED  HOLE 

X=RADIUS*COS ( THETA+ALPHA) 

Y=RAD I US  *  S IN ( THETA+ALPHA ) 

CALCULATE  PARAIIETERS  FOR  LOADED  HOLE  EQUATIONS 

Z1=X+R1*Y 

Z2=X+R2*Y 

Z=X+COMPLX*Y 

MAPPING  FUNCTION 


XXIl=C5QRT(Zl*Zl-DIA*DIA/4.-Rl*Rl*DIA*DIA/4. ) 
XXI2=CSQRT(Z2*Z2-DIA*DIA/4.-R2*R2*DIA*DIA/4. ) 

CHOOSE  CORRECT  SIGN  OF  CSQRT 

90  CONTINUE 
XI1=Z1+XXI1 
XI2=Z2+XXI2 

XI 1=2 . *XI 1 / ( DIA* ( 1 . -COMPLX*R* ) ) 

XI2=2 .*XI2/ ( DIA* ( 1 .-COMPLX*R2) ) 

COXl=REAL(XI 1 ) *REAL(XI1 )+AIMAG(XIl ) *AIMAG (XI 1 ) 
COX  2=REAL ( X 1 2 ) *  REAL ( XI 2 ) +AIMAG ( XI 2 ) * AIMAG ( XI 2 ) 
IF(COXl.GE. .99999)  GO  TO  100 
XXI1=-XXI1 
GO  TO  90 
100  CONTINUE 

IF(COX2.GE. .99999)  GO  TO  110 
XXI2=-XXI2 
GO  TO  90 
110  CONTINUE 
XXI1=XI1 
XXI2=XI2 

CALCULATE  PHI  PRIME 


COM1=(0. ,0. ) 

COM2=(0 . ,0. ) 

DO  120  M=l,45 

COt11»COMl+M*Al  (M)*XI !**(-!  *M) 
COM2=COM2+M*A2 ( M) *XI2** (-1*M) 
120  CONTINUE 
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CHECK  SIGN  OF  CSQ?IT 

XH=CSQRT(Zl*Zl-DIA*DIA/4.-DIA*DIA*Rl*Rl/4.  ) 
XI2=CSQRT(Z2*Z2-niA*DIA/4.-DIA*DIA*R2*R2/4. ) 

CHECK1=Z1/XI1 

CHECK2-Z2/XI2 

IF ( REAL ( CHECK 1) .LT.-. 00001)  XI1=-1.*XI1 
IF(RnAL{CHi:CK2)  .LT.-. 00001)  XI2=-1.*XI2 
PHI1=(AK1-C0H1 )/XIl 
PHI2=(AK2-COM2)/XI2 

CALCULATE  STRESS  COMPONENTS  IN  LAMINATE  AT  COORDINATES  X,Y 

STRX=2 .*REAL{R1*R1*PHI1+R2*R2*PUI2) 

STRy=2.*REAL(PHH+PIII2) 

STRXY=-2 .*REAL(R1*PHI1+R2*PH12) 

STRESS(  1 ,  JJ,NH)=STRX*COS(  ALPH)  *COS(  ALPtl ) +STRY*S IN(  ALPH  )  * 

1  S I N  ( ALPH  )  -  2  .  *STRX  I'*  S I N (  ALPH )  *COS  ( ALPEl ) 

STRESS(2, JJ,NN)=STRX*SIN(ALPH)*SIN{ALPH)+STRY*COS(ALPH)* 

1  COS ( ALPH )+2.*STRXY*SIN{ ALPH )*COS( ALPH) 

STRESSO,  JJ.NN)  =iSTRX*SIN(  ALPH)  *COS(  ALPH) -STRY*SIN(  ALPH)* 

1  COS ( ALPH ) +STRXY* ( COS \ ALPH ) *COS ( ALPH ) - 

2  STN(  Ar,PH)  *SIN(ALPH)  ) 

CALCULATE  DISPLACEMENTS 

XI1=XXI1 
XI2=XXI2 
COM1=(0. ,0. ) 

COM2=(0. ,0. ) 

DO  130  M=l,45 

COMl=COMl+Al (M) *XI1** (-1*M) 

COH2=COM2+A2(M) *XI2**(-1*M) 

130  CONTINUE 

XXI1=CL0G< Xll ) 

XXT.2=CL0G(XI2  ) 

PHI1=AK1*XXI 1+COMl 
PHI2=AK2*XXI2+COil2 
U( JJ,HN)=2.*REAL(P1*PHI1+P2*PHI2) 

V(JJ,NN)  =  2  .*REAL(Q1*PE1I1+Q2*PHI2) 

140  CONTINUE 
150  CONTINUE 

RETURN 

END 
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SUBROUTINE  PLYSTR( IFAIL) 

TRANSFORMS  LAMINATE  STRAINS  TO  PLY  STRESSES/ STRAINS  BY  ASSUMING 
CONSTANT  STRAIN  THROUGH  THE  THICKNESS 

COMMON/ ONE/ El ( 3 ) , E2 ( 3 ) .G12 ( 3 ) , V12 { 3 ) 

COMMON/TWO/10UT(15)  .NUMPLY.NUMMAT,  ANO(8)  .PLYTHKO)  ,MATID(8) 
COMMON/THREE/  lANG,  ILOW  ,  IHIGH,  STPINK,  NUMSTP 
COMMON/ EIGHT/ STRESS ( 3 ,7,73), STRAIN( 3 ,7,73) 
COMMON/NINE/STRl(8,7,73),STR2(8,7,73) ,STR12{8,7,73) 

REAL  IANG,ILOW,IHIGH 


i 


STRAINS  PER  PLY 
MOVE-0 

NUMPT-(  (IHIC,H-ILOW)/IANG)  +  l 
IF(PUTOUT(IOUT, 7) ,EQ.2. )  WRITE(6, 10) 

10  FORMAT(///20X, ‘STRAINS  PER  PLY*,//'  DIST  ANGLE  PLY', 

1'  STRAIN  1  STRAIN  2  SHEAR  STRAIN'/) 

20  CONTINUE 

DO  40  JJ-1, NUMSTP 
DO  40  NN«1,NUMPT 
DO  40  L-1,NUMPLY 
D»ANG(L)*3. 1415926535/180. 0 
STRANX-STRAINd,  JJ,NN) 

STRANY«STRAIN( 2 , J J , NN ) 

GAMA-STRAI N ( 3 , J J , NN ) 

STRANl-STRANX*COS(D)’*COS(D) 

STRAN2-STRANY*  SIN(D)*SIN(D) 

GAMAl 2-GAMA*  SIN(D)*COS(D) 

STRKL,  JJ,NN)-STRAN1+STRAN2+GAMA12 
STRAN1-STRANX*SIN(D)*SIN(D) 

STRAN2-STRANY*COSt  D) *COS( D ) 

GAMA12— 1  .•GAMA*SIN(D)*COS(D) 

STR2 ( L, J J , NN ) - STRAN1+STRAN2+GAMA1 2 
STRANl— 2.*STRANX*SIN(D)*COS(D) 

STRAN2-2 .*STRANY«SIN(D)*COS(D) 

GAMAl 2-GAMA*  COS ( D ) *COS I D ) -GAMA*  SIN(D)*SIN(D) 

STR12(L, JJ,NN)-STRAN1+STRAN2+GAMA12 

ANGLE- (NN-1 )*IANG+IL0W 
DIST-( JJ-1 ) *STPINK 

IP( PUTOUT ( lOUT , 7 ) . EQ . 2 . )  WRITE( 6 , 3 0 ) DIST, ANGLE , ANG ( L ) , 
ISTRKL,  JJ,NN)  ,STR2(L,  JJ,NN),STR12(L,  JJ,NN) 

30  FORMAT(P6.3,2F8.2,3F12.6) 
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STRESSES  PER  PLY 

IF(PUTOUT(IOUT,Q) .EQ.2. )  WRITE(6,50) 

50  FORtdAT(///20X, 'STRESSES  PER  PLY',//*  DIST  ANGLE  PLY', 

1'  STRESS  1  STRESS  2  SHEAR  STRESS'/) 

DO  70  JJ-1,NUMSTP 
DO  70  NN-1,NUMPT 
DO  70  L-1,NUMPLY 
M»MATID(L) 

V21-V12(M)*E2(M)/El(M) 

DEN-1. -V12(M)*V21 
ABC»STR2(L, JJ,MN)/DEN 
BCA-STR1{L, JJ,NN) 

STR1(L, JJ,NM)-E1(M)*STR1(L, JJ,IIN)/DEN+V12(M)*E2(M)*ABC 
STR2(L,  JJ.NH)-V12(.'1)  *E2(M)  ♦BCA/DEH+E2(M)  *A'iC 
STH12(L,  JJ,NN)-STR12(L,  JJ,NIJ)*G12(M) 

ANGLE’* (  Nil- 1  )  *IANG+ILOW 
DlST-( JJ-1)*STPINK 

IF( PUTOUT(IOUT,a) .EQ.2.)  VIRITE( 5, 60)DIST, ANGLE, ANG(L) , 
ISTRKL,  JJ.NN)  ,STR2(L,  JJ,NN)  .STR12(L,  JJ,NN) 

60  FORMAT{F6.3,2Fa.2,3F12.2) 

70  CONTINUE 
MOVE-1 

IFdFAIL.EQ.  1)  GO  TO  20 
80  CONTINUE 

RETURN 

END 


SUBROUTINE  FAILURE 

POINT  STRESS/STRAIN  ANALYSIS  FOR  FAILURE  USING  UNIDIRECTIONAL 
MATERIAL  ALLOWABLES 

COMMON/  TVIO/ 1  OUT  { 1 5 )  ,  NUMPLY ,  N  UMMAT,  ANG  (  8 ) ,  PLYTHK  (  3 )  ,  MATI D  (  8  ) 
COMMON/THREE/ lANG, ILOW, THIGH, STPINK,NUMSTP 
COMMON/ FOUR/ PX, PY, PXY, P, PW, ALPHA, BETA, DIA, CORRECT 
COMMON/ FIVE/ FXT ( 3 ) , FXC ( 3 ) , FYT { 3 ) , FYC ( 3 ) , FXY( 3 ) , I FAI L 
COMMON/NINE/STRl(8,7,73) ,STR2(8,7,73) ,STR12(8,7,73) 

REAL  lANG,  ILOVMHIGIl 

DIMENSION  PLYFAII.(  3,8),  FAIL5{  3 , 0  )  ,  RAT  10 (  3,8),  PLYRATO(  3,8) 

IF(PUTOUT(IOUT,9) .NE.2.)  GO  TO  40 
IF(IFAIL.GT.2)  GO  TO  20 
VIRITE(6, 10) 

10  FORMAT(///2 OX, 'FAILURE  CRITERIA  PER  PLY',//'  DIST  ANGLE  ', 
1'  PLY' , 12X, 'FAILURE  NUMBERS' /35X, • 1  2  SHEAR') 

GO  TO  40 
20  CONTINUE 
WRITE(6, 30) 

30  FORMAT(///20X, 'FAILURE  CRITERIA  PER  PLY',//'  DIST  ANGLE*, 
1*  PLY  FAILURE  FAILURE  RATIOS ' /32X' NUMBER 

2*  1  2  SHEAR') 

40  CONTINUE 


40 


I 


o n n  no  non  on  non  o  on 


CIIECK-0 . 

KKK>-1 
F2-0 . 0 
F3-0.0 

NUMPT-(  (lHIGK-IliOW)/IAK  )♦! 
DO  170  JJ-l.NUMSTP 
DO  170  KK-l.NUMPT 
SIO-1.0 

DO  160  II-l.NUMPLY 

X-PTRKII.JJ.KK) 
Y-STR2(II,JJ,KK) 
XY-STR12(II, JJ.KK) 
MATII-MATIO(II) 


GO  TO  (50.50.60,70,80)  IFAIL 
50  CONTINUE 

MAXIMUM  STRESS/ STRAIN 
FX-PXT(MATII) 

IF(X.LT.O.O)  FX-FXC(MATII) 
FY-PYT(MAIII) 

IF(Y.LT.O.O)  FY-FYC(MATII) 

Fl-X/FX 

F2-Y/FY 

F3-XY/FXY(MATII) 

GO  TO  90 


60  CONTINUE 

TSAI -HILL  CRITERIA 
PX«PXT(MATII) 

IF(X.LT.O.O)  FX-FXC(MATII) 

FY-FYT(MATII) 

IF(Y,LT.0.0)  FY-FYC(MATII) 

F1-X*X/ ( FX*FX)+Y*Y/ (FY*FY)-X*Y/ ( FX*FX)+ 
1  XY*XY/(FXY{MATII)*FXY(^tATII)  ) 
RATIOX-(X/FX)/SQRT(Fl ) 
RATIOY»(Y/FY)/SQRT(Fl ) 

RAT lOXY- ( XY/ FXY ( MATI I ) ) / SQRT( F) ) 

GO  TO  90 


70  CONTINUE 

MODIFIED  TSAI-V/U  CRITERIA 

Fl-1 ./FXT(MATII)-1 ./FXC(MATII) 

F2-1 ./ FYT  ( MATH ) -1  ./ FYC  ( MATH  ) 
F11-1./(FXT(MATH)*FXC(MATH)  ) 

F22-1 ./  (FYT(MATH)*FYC(MATH)  ) 

F66-1 . / { FXY ( MATI I ) • FXY ( MATI I ) ) 

F1*F1*X+F2*Y+F11*X*X+F22*Y*Y+F66''XY*XY 

FX-FXT(MATH) 

IF(X.LT.O.O)  FX-FXClMATII) 
FY«FYT(MATH) 

IF(Y.LT.O.O)  FY-FYC(MATIH 
Fa-ADS(FI) 

RATIOX-( X/FX ) /SQRT( FI ) 

RAT  I OY- ( Y/ FY  )  /  SO  RT  ( F 1 ) 
RATIOXY-(XY/FXY(MATII) )/SQRT(Fl) 

GO  TO  90 
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BO  CONTINUE 


UOFFMAN  FAILURE  CRITERIA 

F1-1./FXT{MATII)-1./FXC(MATII) 

F2«l ./FYT(MATII)-1./FYC(MATII) 

FI  l«l./(FXT(  MATH  )*FXC(  MATH)  ) 

F22-.\  .  /  (  FYT(  MATH )  *FYC  ( MATH  )  ) 

F66-1 . / { FXY( MATI I ) *FXY( MATI I ) ) 

F12— 1  ./(FXT{MATH)*FXC{MATH)  ) 

F1-F1*X+F2*Y+F11*X*X+F22*Y*Y+F12*X*Y+P66*XY«XY 

FX-FXT(MATH) 

IF(X.LT.O.O)  FX-FXC(MATII) 

FV-FYT(MATH) 

IF(Y.LT.O.O)  FY-FYC(MATII) 

Fl-ABS(Fl) 

RATIOX-(X/FX)/SQRT(Fl ) 

RATIOY-(y/FY)/SQRT(Fl) 

RATIOXY- ( XY/ FXY ( MATI I ) ) / SQRT( FI ) 

GO  TO  90 


90  CONTINUE 

ANGLE- (KK-1 ) *IANG+ILOW 
DIST-( JJ-1 ) *STPINK 
IF(IFAIL.GT.2)  GO  TO  110 

IF(PUTOUT(IOUT,9) .EQ. 2. )  WRITE{6,100)  DIST, ANGLE, ANG( H ) , 
IPl «  F2  #  F3 

100  FORflAT(F7.3,2F10.2,3F10.3) 

GO  TO  130 
110  CONTINUE 

IP(PUTOUT(IOUT,9) .EQ.2.)  V;RITE(6, 120)  DIST, ANGLE, ANG( I I ) , 
IFl, RATIOX, RATIOY, RATIOXY 
120  PORMAT(F7.3,2F10.2,4F10.3) 

130  CONTINUE 

AUTOMATIC  SEARCH  FOR  FAILURE  ROUTINE 

IF(Sxa.EQ.2.)  FAII,S(1,H)-F1 

IF(SIG.EQ.2.  )  FAILS(2,  H)-F2 

IF(SIG.EQ.2.  )  FAILS(3,  H)-F3 

IF{SIG.EQ.  2.  )  RATIOU  II)*RATI0X 

IF(SIG.EQ. 2. )  RATIO( 2, H) -RATIOY 

IF(SIG.EQ.2.  )  RAT7:0(  3,  H)-RATIOXY 

IF(JJ.NE.2)  GO  TO  150 

PLYt'AIL(l,H)-Fl 

PLYFAIH  2,  H)-F2 

PLYFAIL(  3,  II)'^F3 

PLYRATO( 1 , I I ) -RATI OX 

PLYRATOC  2 , 1 1 ) -RATIOY 

PLYRATO(  3,  H)-KATIOXY 

CIIK-CHECK 

IF( ABS( CHECK ) .LT. AnS(Fl) )  CHECK-Fl 
IF(An3(CHECK) .LT. AUS(F2) )  CHECK-F2 
IF (AOS (CHECK) .LT. AnS( F3) )  CHECK-F3 
IF { CHECK. EQ.CHK)  GO  TO  150 
KKK-Ki; 

III-II 

SIG-2.0 

DO  140  rl-l,HI 

DO  140  N-1,3 

FAILS(N,  M)=«PLYFAIL(N,  M) 

RATIO(.l,  M)-PLYRATO(N,  M) 

140  CONTINUE 
150  CONTINUE 
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160  CONTINUE 
170  CONTINUE 


IF ( CHECK. RQ. 0.0)  GO  TO  260 

IF(  IFAIL.EQ.  1)  CORRF.CT-1 .0/ABS(CnF.CK) 

IF( IFAIL.EO. 2 )  CORRECT-1 .0/ABS( CHECK) 

IF{  IF/ML.  EQ.  3)  CORRECT-1 . 0/SQRT( CHECK  ) 

IF(IFAIL.CQ.4)  CORRECT-1 . 0/SQRT( CHECK ) 

IFdFML.EQ.S)  CORRECT-1 .0/SQRT(  CHECK) 

IF(PUT0UT(10UT, 10) .HE.2. )  GO  TO  250 
IF(CORRECT.LT. .990. OR. CORRECT. GT. 1 .001 )  GO  TO  250 
ANGLE- ( :<KK-1  )*IANG+ILOv; 

IF(1FAIL.GT.2)  GO  TO  210 
WRITfUG,  1^0)  PX,PY,PXY,P 

ino  FORtlAT(///20X, 'AUTOMATIC  SEARCH  BOR  FAILURE* //25X, 

1' FAILURE  STRESSES' / IBX, ' PX' , lOX, ' PY' , lOX, ' PXY* , lOX, ' P' /, 
211X,4F12.2//9X, 'DIST  ATJGLE  PL^i  n  FAILURE  NUMBERS' / , 

336X,'l  2  SHEAR') 

DO  200  I-1,NUMPLY 

WUITEC  6 , 190 )  STPI NK, ANGLE, AHG( I ) , FAILS (1,1), FAILS( 2 , 1 ) , FAILS ( 3,1) 
190  PORMAT(Fll.3,P9.2,F10.2,3F0.3) 

200  CONTINUE 
GO  TO  260 

210  CONTINUE 

NRITE(6,220)  PX,PY,PXY,P 

220  POR;1AT(///20X, 'AUTOMATIC  SEARCH  FOR  FAI  LURE"  /  /  2  5X, 

1 'FAILURE  STRESSES' /ISX, ' PX' , lOX, ' PY* , lOX, ' PXY' , lOX, ' P' / , 
211X,4F12.2//5X, '  DIST  ANGLE  PLY  FAILURE  ', 

39X, ' FAILURE  RATIOS' /35X, 'NUMBER' , 7X , ' 1 ' , 9X , ‘ 2 ' , 7X , ' SHEAR' ) 

DO  240  I-1,NUMPLY 

V;RITE(&, 230)  STP INK, ANGLE, AMG( I ) , FAILS! 1 , I ) , RATIO ( 1,1), 

1  RATIO(2,I) , RATIO! 3,1) 

230  FORMAT(F11.3,F9. 2,F10.2,4F10.3) 

240  CONTINUE 
250  CONTINUE 


260  CONTINUE 

IF! PUTOUT! IOUT, 10) .NE.2. )  CORRECT-1. 0 

RSrURN 

END 
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sun rtoirr  i  n  t-:  s  imult  (  a  ,  u  ,  •  i  ,  ks  ) 

T:i:nT  ?o(!;  AiiOoniTiMic  s:,':  fGUL.\Rmr  »,om:D  oi/io/T9 
DIMEMSION  A(1).U{1) 

MACiii:  rj  ':;paiirjON  for  cyijer  stwee  rRECisioti 
DATA  GP':i/7.1LE-15^ 
roL-:t.  i:5*r,PS'*(M-i ) 

RETA-'O .  O 

J  J*- ?I 

00  r.s  j»i,N 

JY<T+1 
0J-JJ+-i  +  l 
ijT(',A«n.o 
IT-JJ-J 
lOO  30 
IJ-IT+I 

IF(AM{  r,I-3A)-AriS(  A(  IJ)  )  )  20,30,30 
20  ,1  IGA- A (1.7) 

IMAX-I 

30  co>7Ti  iu:; 

IF(  Ans(  3IGA)  .GT.  nF.TA);3ET.\-A3S(QIGA) 
lF(AnS(  OIGA)-T()L*DnTA)  3S,  35,40 
35  t;S-1 
RETURN 

40  H-J+:l*(J-2) 

IT-IMAX-J 
L'O  50  K»J,N 

11- Il+N 

12- Il+IT 
SAVl>A(  1 1 ) 

A(ti)-A(I2) 

A( I2)-SAVE 

50  A(Il)«A(Il)/aiGA 
oAVl>3(  I  MAX) 

U(i:iAX)-0(  J) 

B( J)-SAVE/BIGA 
IV(J-N)  55,70,55 
55  I0S-N»(J-1) 

DO  65  IX-JY,N 

IXJ-IQS+IX 

IT-J-IX 

DO  GO  JX»JY,N 

IXJX-N*(JX-l)+IX 

j  j;<»ixJx+iT 

60  A(  IXJX)-A(  IX.7X)-{A(  IXJ)*A(JJX)  ) 

65  0(IX)»B(IX)-(3(J)*A(IXJ)) 

70  JY-N-l 

DO  00  J»l,iIY 

];a=iIT-j 

IC»N 

DO  90  K>1,J 

a(  in)-ii(iD)-A(  iA)*o(  ic) 
lA-IA-M 
80  IC-lC-l 
RETURN 
END 
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SUnROUTlNE  ROOTS (XCOF.COF.M.ROOTR, ROOTI, IKR) 
DIHKNSION'  XCOF(M)  ,C0F(M)  ,  R0(7TR(  M)  ,  ROOTI  (  M) 

DOUBLE  PRECISION  XO, YO, X, Y, XPR. YPR . UX. UY, V , YT , XT , U 

1  ,XT2,YT2,SU'‘iSO,  DX,DY,TEMP,  ALPHA,  FI 

2  ,RMPRF,C,TOL 

C  RELATIVE  MACHINE  PRECISION  (TEST  FOR  'ALMOST  ZERO’) 

DATA  RMPREC/l .00-14/  ,TOL/l.OD-4/ 

IFIT-0 

N-M 

IER-0 

IF(XCOF(N+l ) )  10,30,10 
10  IF(N)  20,20,50 
20  IER-1 

GO  TO  360 
30  i::r-4 

GO  TO  360 
40  IER-2 

GO  TO  360 

50  IF{N-36)  60,60,40 
60  NX-ri 
NXX-M+1 
N2-1 

KJl  »  N41 
DO  70  L-1,KJ1 
MT-KJl-L+l 
70  COF(MT)«XCOF(L) 

90  XO-.00500101D0 
Yo-o.oioooioino 

IN-0 
90  X-XO 

XO— IO.DO*YO 
YO— 10.D0*X 
X-XO 
Y=YO 
iN-i:i+i 
GO  TO  110 
100  IFIT-1 
xrR».x 

YPR=-Y 
110  ICT-0 
120  UX-O.DO 
UY»0 . 00 
V  -0  .  DO 
YT-o.no 
XT=1 .00 
U=Cf)F(N  +  l  ) 

IF  (  DAHS  ( U  )  .  Ur. .  RMPREC )  GO  TO  200 
130  00  140  1=1, N 
L  =;i-i  +  i 
TUMP-COF(L) 
y.T2-X*XT-Y*YT 
YT2-X*YT+Y*XT 
U-U+TEMP*XT2 
V-V+TEMP*YT2 
FI -I 

UX=UX+FI'*XT*TEMP 
UY-UY-FI*YT*TRMr 
XT-XT2 
140  YT-YT2 

SUMSO-UX*UX+UY*UY 
1F(SHMSQ.LE. RMPREC)GO  TO  230 
150  DX-(V*UY-U*UX)/SUMSQ 
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X-X+DX 

DY—  ( U*UY+V*UX  )  / SUMSQ 
Y-Y+DY 

160  IF(  DAnS(DY)+  DADS( OX) .LE.TOL)  GO  TO  210 
170  ICT-ICT+1 

IF(ICT-500)  120,180,ia0 
180  IF(IFIT)  210,190,210 
190  IF(IN-5)  90,200,200 
200  IRR-3 

GO  TO  360 

210  DO  220  L»1,NXX 
;1T-KJ1-L+1 
RTRMP-XCOFiMT) 

XCOF(MT)-COF(L) 

220  COF(L)-RTEMP 
ITEMP-N 
;j*MX 

NX-ITEMP 

IF(IFIT)  250,100,250 
230  IP(IFIT)  240,90,240 
240  X-XPR 
V»YPR 
GO  TO  210 
250  IFIT-0 

260  IF(DAn3(Y)-1.0D-4*DABS(X) )  290,270,270 
270  ALPMA*X+X 

SUMSQ-X*X+Y*Y 
N»N-2 
GO  TO  300 
290  X»0.00 
tIX«MX-l 
NXX-rr<X-l 
290  Y-O.DO 

SUMSQ-O.DO 

ALPHA=*X 

N-tI-1 

300  COF(2)=COF(2)+ALPHA*COF(1) 

310  DO  320  L»2,N 

320  COF ( L+1 ) -COF ( L+1 ) +ALPHA*COF ( L) -SUMSQ*COF ( L- 1 ) 
330  ROOTI(N2)-Y 
ROOTR(N2 )-X 
N2-M2+1 

IF(SUMSQ.I  R. RMPRCC)  GO  TO  350 
340  Y»-Y 

SUnSQ»0.D0 
GO  TO  330 

350  IF(N.GT.O)GO  TO  80 
360  RETURN 
END 
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